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ABSTRACT 

We make a detailed investigation of the properties of Lyman-break galaxies (LBGs) in the 
ACDM model. We p resent predictions f or two published variants of the GALFORM semi- 
analytical model: the iBaugh et al.l (120051) model, which has star format ion at high redshifts 
dominated by merger-driven starbursts with a top-heavy IMF, and the iBower et al.l ([2006) 
model, which has AGN feedback and a standard Solar neighbourhood IMF throughout. We 
show predictions for the evolution of the rest-frame far-UV luminosity function in the red- 
shift range z = 3 — 20, a nd compare with the observed luminosity functions of LBGs at 
z = 3 — 10. W e find that the Baug h et all model is in excellent agreement with these observa- 
tions, while the lBower et al.l model predicts too many high-luminosity LBGs. Dust extinction, 
which is predicted self-consistently based on galaxy gas contents, metallicities and sizes, is 
found to have a large effect on LBG luminosities. We compare predictions for the size evo- 
lution of L BGs at different luminosities with observational data for 2 < z < 7, and find the 
Baug h et"ai] model to be in good agreement. We present predictions for stellar, halo and gas 
masses, star formation rates, circular velocities, bulge-to-disk ratios, gas and stellar metal- 
licities and clustering bias, as functions of far-UV luminosity and redshift. We find broad 
consistency with current observational constraints. We then present predictions for the abun- 
dance and angular sizes of LBGs out to very high redshift (z ^ 20), finding that planned 
deep surveys with JWST should detect objects out to z < 15. We predict that the effects of 
dust extinction on the far-UV luminosity density should be large (~ 2 mag), even out to high 
redshifts. The typical UV luminosities of galaxies are predicted to be very low at high red- 
shifts, which has implications for detecting the galaxies responsible for reionizing the IGM; 
for example, at z — 10, 50% of the ionizing photons are expected to be produced by galaxies 
fainter than M AB ( 1500 A) - 5 log h 15. 
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1 INTRODUCTION 



The d iscovery of Lyman-break galaxies at z ~ 3 by Stei del et al.l 
dl996h was a breakthrough in observational studies of galaxy for- 
mation. It provided, for the first time, a significant sample of normal 
galaxies at high redshift whose properties and population statistics 
could then be investigated obs ervationally and com pared to the pre- 
dictions of theoretical models teaugh et al.l lT998). 

Lyman-break galaxies (LBGs) are star-forming galaxies 
which are identified through the Lyman-break feature in their 
spectra. This feature is produced by absorption by neutral hy- 
drogen in the atmospheres of massive stars, in the interstellar 
medium (ISM) of the galax y and in the intergalactic medium (IGM) 
dSteidel & HamiltorJI 19921) . For ground-based observations, detec- 
tion of the Lyman break is restric ted to redshifts z > 3. Since the 
first successful demonstration by ISteidel et alj {T99ft) at z ~ 3, 
the technique has been extended to identify galaxies at both higher 
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redshifts and lower luminosities, using ground -based telescopes 
and the Hubble Sp ace Telescope (HST) (e.g. iMadau et all 1 1996 



Steidel et al] 



Yoshida et al 



19991: iBouwens et alj [20031 : Shimasak uetaiTl200; : 
2006). Finding LBGs at z > 6 requires observing 



in the near-IR, which was first done using NICMOS on HST, and 
led to the identification of a small number of z ~ 7 — 8 objects 
jBouwens et alj |2004b). More recently, the WFC3/IR camera on 
HST has all owed the discovery o f much larger samples of LBGs at 
z ~ 7 - 8 jBouwens et al|20ld:lMcLure et alj|20ld : l0esch etail 
l2010b! : lBunker et al.ll2009ljYan et alj|2009h and even a few candi- 
dates at z ~ 10 (Bo uwens et al. 2009b). By observing in the UV 
from s pace, the GALEX sate llite has also been used to find LBGs at 
z ~ 1 jBurgarella etai]|200rj) . 

Follow-up observational investigations on LBGs have in- 
cluded estimates of their star formation rates (SFRs), sizes, mor- 
phologies, stellar and dynamical masses, galactic outflows, metal- 
licities, dust extinctions, gas masses, IR/sub-mm emission from 



dust and clustering (Steidel et alj 1 19961; iGiavalisco et al.l [T996 
lAdelberger et alj 1 19981 : ICh apman et alj l200d : IPettini et alj 12001 
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Shaple v et alj200ll ; lFerguson et alj2004h . Since LBGs are selected 
on the basis of their rest-frame far-UV emission, which is domi- 
nated by massive young stars, LBG samples at different redshifts 
also provide a means to trace the cosmic star forma tion history, 
a key component in our picture of galaxy formation CMa dau et al.l 
1 19961) . although important uncertainties remain due to the effects of 
dust extinction. Observations of LBGs at z — 3 — 10 probe galaxy 
evolution over the first 3-15% of the age of the universe. 

Other observational techniques have also been used to find 
normal galaxies at high redshift. S earches for star -forming galax- 
ies via their Lya emission line (e.g. lHuetalJll998l) cover a similar 
redshift range to LBGs. The main drawback with this technique 
is that so me star forming-gal axies show Lya absorption rather than 
emission dShaplev et al.l2003l) . and consequently are missed in nar- 
rowband surveys. A further complication is that inferring the SFR 
from the Lya emission line is much more uncertain than inferring 
it from the far-UV continuum, since the effects of dust extinction 
are amplified by resonant scattering of Lya photons by hydrogen 
atoms. Another technique is to search for sub -mm or IR emis- 
sion from dust in h igh-z star-forming galaxies dSmail et aljfl997l ; 
Hugh es et al. I ll998h . This method is currently limited by source 
confusion at faint fluxes due to the relatively poor angular reso- 
lution of current IR/sub-mm telescopes, which restricts searches to 
redshifts z < 3 and mostly to the galaxies with the highest SFRs. 
Other techniques for selecting high-z galaxies, which are sensitive 
also to non-star-forming galaxies (such as ERO and DRG colour 
selection) are limited to even lower redshifts z < 2. Overall, the 
Lyman-break technique still seems the most effective for finding 
large samples of star-forming galaxies at z > 3 that cover a wide 
range of luminosities and SFRs. 

The theoretical significance of the discovery of LBGs was 
highli ghted early on us ing semi-analytical models of galaxy forma- 
tion. iBaughet all d 19981) showed that the abundance and observed 
properties of the z ~ 3 - 4 LBGs found by Steid el et ail d 19961) and 
Mada u et alj dl996t) could be explained in the framework of CDM, 
and that they fitted into a picture in which the cosmic SFR density 
peaked around a redshift z ~ 2. This evolution of the cosmic SFR 
density was driven by the combined effects of the build-up of dark 
matter halos, gas cooling and supernova feedback. Their model had 
star formation occuring mostly in quiescent disks, and neglected the 
effects of dust extinction. Subsequent observational studies found 
evidence from UV continuum slopes for signific a nt dust extinction 
in LBG s dMeurer et alJI 19991: ISteidel et alj]l999l) . ISomerville et al.l 
d200ll) proposed a different semi-analytical model in which star 
formation bursts triggered by galaxy mergers played an important 
role, and combined this with an empirical prescription for dust ex- 
tinction, tuned to match observational estimates of the extinction 
for z ~ 3 LBGs. 

More recent studies of L BGs in semi-analytic al models in- 
clude |Guo|& Whitel d2009t) andlLo" Faro et al. (2009), which inves- 
tigated the effects on inferred luminosity functions and other prop- 
erties of applying observational LBG colour selections to model 
galaxies. The former used a phenomonelogical model for dust ex- 
tinction, while the latter used a physical model similar to that in the 
present paper. LBGs wer e also studied in gas-dynamical simula - 
tions of galaxy formation dNagaminel l2002 : Wei nberg et al.l l2002). 
but these simulations had the drawback that they did not predict 
properties for the present-day galaxy population consistent with ob- 
servations, unlike the semi-analytical models. Furthermore, none 
of these models were able to explain the number counts and red- 
shifts of faint sub-mm galaxies discovered in surveys at 850/im 



dSmailetal.|[l997l) , wh ich were subsequently shown to be dusty 
- 3 Jctiapman et alj2003h . 



starbursts at z 



In order to explain within a single framework the sub-mm 
and Lyman-break galaxies at high redshift, together with a wide 
range of galaxy properties at z = (including optical and near- 
and far-IR l uminosity functions , gas fractions, metallicities and 
galaxy sizes), lBaughetalJd2005l) introduced a new semi-analytical 
model in which the gas consumption timescale in disks at high 
redshifts was increased, with the result that starbursts triggered 
by galaxy mergers played a more significant role at high redshift. 
They assumed, further, that stars formed in these bursts with a 
top-heavy initial mass function (IMF). Unlike previous models of 
high-redshift galaxies, this model included a fully self-consistent 
treatment of both absorption and emission of radiation by dust, 
with dust extinction calculated from radiative transfer based on the 
predicted gas masses, metallicities and sizes of galaxies, and the 
spectrum of the dust emission calculated by solving for the tem- 
perature distribution of the dust grains in each galaxy. In subse- 
quent papers, we have explored other predictio ns from the same 
model , including stellar and gas metallicities (Nagashima et al. 
2005 a.b), g alaxy colours, si z es and mo rphologies in the local uni- 
verse dAlmeida et al.l 120071: iGon zalez et al. 200 9J), the evolutio n 
of Lva-emitters dLe Delliou et alj|2005l, 120061 : lOrsi etal] |2008l) . 
and the evolution of ga laxies at mid- and far-IR wavelengths 
dLacev et alJl2008Ll2010bl) , finding generally good agreement with 
observational data. 



In iBaugh et al.l d2005l) . we made only a limited comparison 
with observational data on LBGs, focussing on their rest-frame 
far-UV luminosity function at z = 3. Since then, there has been 
a huge increase in the amount and quality of observational data 
on LBGs, in particular enabling measurements of the luminosity 
function of LBGs out to z ~ 10. Therefore in this paper we re- 
turn to studying LBGs, making detailed predictions for the evo- 
lution of their luminosity functions over a wide redshift range 
(z = 3 — 20) and for many other properties. We consider tw o 
variants of the GALFORM s emi- analytical model \Cole et alj 2000), 
those of lBaugh et ail d2005h and lBower etal] d2006l) . The two mod- 
els differ in a number of significant ways, the most important be- 
ing that the Bower et"al] model includes AGN feedback, while the 



IBaugh et alj model has a variable IMF, as already mentioned. We 
focus here on far-UV-selected galaxies in the redshift range z > 3, 
where they are observationally detected using their Lyman-break 
features. We inve stigate the present-day descendants of LBGs in a 
companion paper dGonzalez et al]2010al) . and make predictions for 
the reio nization of the IGM from the same models in lRaicevic et al.1 
( 2010a. b). The properties and evolution of far-UV-selected galaxies 
at lower redshifts will be the topic of a separate paper. 

The plan of this paper is as follows: In Sj2] we briefly review 
the main features of the two models. In ij3] we compare predictions 
for the evolution of the far-UV luminosity function with observa- 
tional data from LBGs, and investigate the sensitivity of these pre- 
dictions to various model parameters. In fj4] we investigate the sizes 
and other physical properties of UV-selected galaxies, and carry out 
a detailed comparison with the observed sizes of LBGs. In ij5] we 
present predictions for LBGs at very high redshifts, which may be 
accessible with future telescopes such as JVKSTand ELTs. In ij6] we 
show how our predictions for LBGs fit into the wider picture of the 
evolution of the cosmic star formation and UV luminosity densities. 
We briefly consider the contribution of LBGs to the reionization of 
the IGM. Finally, we present our conclusions in fjT] 
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2 THE GALFORM GALAXY FORMATION MODEL 

We compute the formation and evolution of galaxies within the 
framework of the ACDM model of structure formation using the 
semi-analytical galaxy formation model GALFORM. The general 
methodolog y and approxim ations behind GALFOR M are set out i n 
detail in ICole et ail I l2000h (see also the review by lBaugb]|2006h . 
In summary, GALFORM follows the main processes which shape 
the formation and evolution of galaxies. These include: (i) the col- 
lapse and merging of dark matter halos; (ii) the shock-heating and 
radiative cooling of gas inside dark halos, leading to the forma- 
tion of galaxy disks; (iii) quiescent star formation in galaxy disks; 
(iv) feedback from supernova explosions, from AGN heating, and 
from photoionization of the IGM; (v) chemical enrichment of the 
stars and gas; (vi) galaxy mergers driven by dynamical friction 
within common dark matter halos, leading to the formation of stel- 
lar spheroids, and also triggering bursts of star formation. The end 
product of the calculations is a prediction of the numbers and prop- 
erties of galaxies that reside within dark matter haloes of different 
masses. The model predicts the stellar and cold gas masses of the 
galaxies, along with their star formation and merger histories, their 
disk and bulge sizes, and their metallicities. The stellar luminosi- 
ties of the galaxies are then computed from their star formation and 
chemical enrichment histories using a stellar population model. Fi- 
nally, the dust extinctions at different wavelengths are calculated 
self-consistently from their gas and metal contents and sizes using 
a radiative transfer model. 



2.1 Baugh05 and Bower06 models 

The two variants of GALFORM c onside red in this paper, 
iBaugh et al.l J2005l) and lBower etal.1 J2006l) (hereafter Baugh05 
and Bower06), have been used as fiducial models in a number of in- 
vestigations. We now briefly compare the m ain features of the two 
models. Similar discussions can be found in Almeida et al. (2007, 
2008), Gonzalez et al]j2009l) andlGonzalez-Perez et alj<2009h. For 



full details o f the two mod els, see the p a pers bv lCole et al.l fc OOO). 
IBaugh et al.l d2005h (also lLacev et all J2008h \ and iBower et all 



(2006). We give a summary of previously published results ob- 
tained with the two models, and of their relative performances in 
explaining different observed properties of galaxies, in the Ap- 
pendix. 

(i) Cosmology: The Baugh05 model adopts a ACDM cosmol- 
ogy with a present-day matter density parameter, fl m = 0.3, a 
cosmological constant, Qa = 0.7, a baryon density, Qt> = 0.04, 
a Hubble constant h — 0.7 in units of 100km s~ 1 Mpc _1 , and a 
power spectrum normalization given by og — 0.93. The Bower06 
model instead uses th e cosmological mod el assumed in the Mil- 
lennium simulation (Springel et al] 120051) . where Q m = 0.25, 
Q A = 0.75, n b = 0.045, h = 0.73 and a s = 0.9. 

(ii) Star formation in disks: In both models, stars can form either 
quiescently in disks or in starbursts at a rate ip = M gas /r», where 
M gas is the mass of cold gas and the timescale, r„, is different in 
disks and bursts. The two models adopt different parameterizations 
for the dependence of the quiescent timescale on disk properties. In 
the Baugh05 model, this timescale varies simply as a power of the 
disk circular velocity, while in the Bower06 model it is proportional 
also to the disk dynamical time. Since the typical dynamical time 
gets shorter with increasing redshift, the star formation timescale 
in the Bower06 model is shorter at high redshift than it would be 
in the equivalent disk in the Baugh05 model. As a consequence, 
galactic disks at high redshift tend to be gas poor in the Bower06 



model, but gas rich in the Baugh05 model. This then results in much 
more gas being available to fuel starbursts and in a higher fraction 
of star formation occuring in bursts at high redshift in the Baugh05 
compared to the Bower06 model. 

(iii) Starbursts: In the Baugh05 model, starbursts are triggered 
only by galaxy mergers (both major and minor), while in the 
Bower06 model, starbursts are triggered also by disk instabilities. 
(We define major mergers as those in which the mass ratio of the 
smaller to larger galaxy exceeds 0.3.) Both major mergers and 
disk instabilities are assumed to transform the stellar disk(s) into 
a spheroid, while in minor mergers only the stars from the smaller 
galaxy are added to the spheroid, leaving the stellar disk of the 
larger galaxy intact. Stars formed in bursts are always added to the 
spheroid. The star formation timescale in bursts is assumed to be 



7*burst rri'lX^dyn'Tdyn , T^burst ,n 



(1) 



so that it scales with the spheroid dynamical time Td yn , but with 
a "floor" value when this is very short. The parameters /d yn 
and r.burst.min are different in the two models: /d yll = 50 and 
T*burst,min = 0.2Gyr in the Baugh05 model and /d yn = 2 and 
T"*burst,min = 0.005Gyr in the Bower06 model, so that the burst 
SFR timescale is 25-40 times larger in the Baugh05 model com- 
pared to the equivalent galaxy in the Bower06 model. The SFR in 
a burst then decays exponentially with time, with a timescale that 
depends on the supernova feedback efficiency and gas recycling 
fraction as well as the SFR timescale. The burst is assumed to be 
truncated after n T e-folding times, at whic h point the remainin 
gas in the burst is ejected into the halo (see lGranato et alj d20Q( 
for more details). The Baugh05 and Bower06 models both assume 
n T — 3. 

(iv) Stellar init i al ma ss function (IMF): The Bower06 model 
uses the lKennicutJ ( ll983l) IMF, consistent with deductions from the 
Solar neighbourhood, in all modes of star formation. This has the 
form diV/dlnm oc m~ x (where N is the number of stars and m 
is the stellar mass) with slope x = 0.4 for m < Mq and x = 1.5 
for m > Mq (compared to x — 1.35 for the Salpeter IMF). The 
Baugh05 model also adopts this IMF in quiescent star formation in 
galactic disks. However, in starbursts triggered by galaxy mergers, 
a top-heavy IMF is assumed, with slope x = 0. The IMF covers 
the stellar mass range 0.15 < m < 12OA/0 in all cases. The yield 
of metals (p) and the fraction of gas recycled per unit mass of stars 
formed (R) are chosen to be consistent with the form of the IMF. 
We assume p = 0.023 and R = 0.41 for the Kennicutt IMF, and 
p = 0.15 and R = 0.91 for the top-heavy x = IMF. 

(v) Supernova (SN) feedback: In both models, supernova explo- 
sions are assumed to reheat cold gas (in both disks and starbursts) 
and eject it into the halo, at a rate 



M eject = (V hot /V m l) a ^lP, 



(2) 



where Vgai is the circular velocity of the disk (for quiescent star for- 
mation) or spheroid (for bursts), and Vhot an d Qhot are parameters. 
The SN feedback is much stronger in the Bower06 model (Vhot = 
485kms _1 and <»hot = 3.2, compared with Vhot = 300kms _ 
and ahot = 2 in the Baugh05 model). A further difference is that 
in the Bower06 model ejected gas is reincorporated into the halo on 
a shorter timescale than in the Baugh05 model. This SN feedback 
supresses star formation much more effectively in low mass than 
high mass galaxies. 

(vi) Superwind vs AGN feedback: A major difference between 
the two models is in the modelling of feedback in massive galaxies. 
If only the "standard" SN feedback described above is included, 
then too much gas cools in massive halos, resulting in the break 



© 0000 RAS, MNRAS 000, 000-000 



4 Lacey et al. 



at the bright end of the present-day optical and near-IR galaxy lu- 
minosity functions being at too high luminosity. In the Baugh05 
model, this pr oblem is solved by in trod ucing supernova-dr iven su- 
perwinds (see lBenson et al] d2003h and lLacev et al.l {2008) for de- 
tails), which eject gas completely from the halo at a rate propor- 
tional to the SFR, and operate alongside the standard SN feed- 
back. Such winds have been observed in massive galaxies, with 
the inferred mass ej ection rates found to be comparable to th e star 
formation rate (e.g. iHeckman et alJI 19901 : IPettini et ai1l200lh . The 
effect of expelling gas from halos is to reduce the density of hot 
gas in them and increase the radiative cooling time, resulting, in 
particular, in less gas cooling in massive halos. In contrast, the 
Bower06 model solves the same problem by introducing feedback 
from AGN. Halos in which the cooling time of the gas exceeds 
the free-fall time and which also contain a sufficiently large central 
supermassive black hole (SMBH) are assumed to set up a steady 
state in which halo gas accreting onto the SMBH releases energy 
in the form of relativistic jets which heat the halo gas, exactly bal- 
ancing the halo cooling, and preventing any gas from cooling onto 
the galaxy disk. The SMBH are assumed to grow by gas accretion 
in starbursts triggered by galaxy mergers and disk instabilities and, 
less commonly, by black hole/black h ole mergers (as detailed i n 
iBower et alj|200d iMalbon et alj|2007l and iFanidakis et alj|2010h . 
In the Bower06 model, the disk instabilities appear to be critical 
for producing large enough SMBH at early enough times for the 
AGN feedback to be effective. 

The assumption of a top-heavy IMF in starbursts in the 
B augh05 mode l is a controversial one. As discussed in detail 
in Bau gh et al.l d2005h . this top-heavy IMF was found to be re- 
quired in order to reproduce the observed number counts and red- 
shift distributions of the faint sub-mm galaxies. The top-heavy 
IMF results both in higher bolometric and far-UV luminosities 
for young stellar populations, and greater production of heavy ele- 
ments and hence also dust, both effects being important for repro- 
du cing the properties of SMGs in the model. Furthermore, as shown 
bv lNagashima et al] d2005al lbl). the predicted chemical abundances 
of the X-ray emitting gas in galaxy clusters and of the stars in ellip- 
tical galaxies also agree better with observational data in a model 
with the top-heavy IMF in bursts, rather than one with a universal 
Solar neighbourhood IMF. Subsequent work using the same model 
has also shown that it predicts galaxy evolutio n in the mid-IR i n 
good agreement with observations by Spitzer dLacev et alj|2008l) . 
A more detailed comparison of the mod el with the properties o f 
observed SMGs has been carried out by ISwinbank et alj d2008t) . 
and an investigation o f the present-d ay descendants of SM Gs in 
Gonzalez etalTd2010bh . As shown bv lLe Delliou et alj d2006h and 
Orsi et al.l (2008), the same model also reproduces the observed 
evolution of the luminosity function and clustering of Ly«-emitting 
galaxies at high redshift. 

A variety of other observational evidence has accumulated 
which suggests that the IMF in some environments may be top- 
heavy compared to the Sol ar neighbour h ood IM F (see lElmegreenl 
2009 for a recent review). Rie ke et al] dl993l) argued for a top- 
heavy IMF in the nearby s tarburst M82, bas ed on modelling its 
integrated properties, while [parra et al] d2007l) found possible evi- 
dence for a top-heavy IMF in the ultra-luminous starburst Arp220 
from the relative numbers of supernovae of different types observed 
at radio wavelengths. Evidence has been found for a top-heavy IMF 
in so me star clusters in int ensely star-forming regions, both in M82 
(e.g.lMcCradv et al.l2003l). and in our own Galax y (e.g. lFiger et alj 
1 19991 : IStolte et al.l 120051 : lHaravama et al]|20o3) . Observations of 



both the old and young stellar populations in the central 1 pc of 
our Galaxy also f avour a top-heavy IMF (Paumar d et al] 1 20061 : 
iManess et alj2007l) . In the local Universe, Meur er et al] d2009l) find 
evidence for variations in the IM F between galaxies from variations 
in the Ha/UV luminosity ratio. Fardal et al. (2007) found that rec- 
onciling measurements of the optical and IR extragalactic back- 
ground with measurements of the cosmic star formation history 
also s eemed to require an averag e IMF that was somewhat top- 
heavy. |Pere z-Gonzale z~et alj d2008l) compared observational con- 
straints on the evolution of the star formation rate density and stel- 
lar mass density over cosmic time, and found that reconciling these 
two types of data also favours a more top-heavy IMF a t higher red- 
shifts, as had been hinted at by earlier studies. Finallv. lvan Do kkum 
d2008l) found that reconciling the colour and luminosity evolution 
of early-type galaxies in clusters also favoured a top-heavy IMF. 
ILarsorJ dl998l) summarized other evidence for a top-heavy IMF 
during the earlier phases of galaxy evolution, and argued that this 
could be a natural consequence of the temperature depend ence of 
the Je ans mass for gravitational instability in gas clouds. ILarsorJ 
d2005h extended this to argue that a top-heavy IMF might also be 
expected in starburst regions, where there is strong heating of the 
dust by the young stars. 

2.2 Dust extinction. 

The attenuation of starlight by dust is a crucial part of our model, 
since LBGs are detected by their rest-frame far-UV emission, 
which appears to be heavily e x tincte d. In our earlier investi- 
gation of LBGs in iBaugh et alj d20050 . we calculated both the 
absorption and emission of starlight by dust in model galaxies 
by coupling GAL FORM to the GRASIL spectrophotometric code 
(Silvaetal. 1998), which incorporates a radiative transfer calcu- 
lation of starlight through a realistic dust distribution, and a de- 
tailed model of heating and cooling of dust grains (including PAH 
molecules). The GRASIL model assumes a composite disk plus 
spheroid geometry for quiescent galaxies, with stars in both com- 
ponents, but dust only in the disk. In the case of ongoing bursts, 
this is replaced by a stellar spheroid plus a flattened burst com- 
ponent containing the dust, and having the same half-mass radius 
as the spheroid. The dust is assumed to be in two phases, a dif- 
fuse medium and also giant molecular clouds in which the stars 
form, and from which they escape after a few Myr. The coup led 
GALFORM+GRASIL model is described in lGranatoetai]d2000l) . 

The GRASIL model is computationally expensive, since it 
involves a detailed radiative transfer calculation for each model 
galaxy, and furthermore in the present paper we only need to calcu- 
late the extinction of starlight from galaxies, and not the IR/sub-mm 
dust emission. In this paper, we therefore calculate dust extinction 
using an alternative approach, which assumes the same geometry 
for the stars and dust as in GRASIL, but gains in computational 
speed by making some simplifying a pproximations. T his alterna- 
tive method is described in detail in lLacev et alj d2010al) . where we 
also show that it gives very similar results to GRASIL for the dust 
extinction and for the far-IR (but not mid-IR) emission from dust. 
Here we give a brief summary of how dust extinction is calculated 
in this alternative approach. 

The starting point of the approach used here is the model 
for extinct i on by diffuse, smoothly distributed dust described in 
Cole et al. (2000), to which we add a treament of extinction of 
young stars by molecular clouds. We calculate the extinction by the 
diffuse compo nent of the dust usin g the tabulated radiative trans- 
fer models of iFerrara et alj dl999l) . which assume that the stars 
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are in both a disk and a spheroid, and that the dust is uniformly 
mixed with the stars in the exponential disk, with a Milky Way 
extinction law. This already gives very different extinctions from 
a fore ground scre en model with the same column density of dust. 
The iFerrara et alj models give dust attenuation factors separately 
for the disk and spheroid light, as functions of wavelength, galaxy 
inclination, central face-on V-band dust optical depth, and ratio of 
spheroid to disk scalelengths. 

The radii of the disk and spheroid of each model galaxy ar e 
predicted directly by GALFORM, as described in lCole et ail {2000), 
based on conservation of angular momentum (for the disk) and en- 
ergy (for the spheroid), and ta king account o f their own self-gravity 
and that of the dark halo. The lFerrara et al.l models assume a fixed 
ratio of vertical to radial exponential scalelengths for the disk stars 
of 0.1. We assume that the dust has the same radial and vertical 
scalelengths as the stars. The total mass of dust in each galaxy 
is calculated from the cold gas mass and metallicity predicted by 
GALFORM, assuming that the dust-to-gas ratio is proportional to 
the gas metallicity, and normalized to match the solar neighbour- 
hood dust-to-gas r atio for solar met allicity (Zq = 0.02). The only 
change relative to lCole et alj d2000h is that we assume that a frac- 
tion /doud of the cold gas and dust are in molecular clouds, so the 
mass of dust in the diffuse phase is a fraction 1 — / c ioud of the total 
dust mass. We then calculate the central face-on optical depth of the 
diffuse dust from its mass and radial scalelength. Finally, we choose 
a random inclination angle for each model galaxy. GALFORM cal- 
culates separate disk and spheroid luminosities at each wavelength, 
so we interpolate in the tables and then apply the separate dust at- 
tenuation factors for the disk and spheroid light. 

We model extinction by molecular clouds based on the same 
physical picture as in GRASIL, but by a more approximate tech- 
nique. We assume that a fraction / c i ou d of the cold gas and dust is 
in molecular clouds of mass M c i ou d and radius r c ioud- All stars are 
assumed to form inside clouds, and then to escape on a timescale 
tesc, such that the fraction inside clouds drops continuously from 1 
for ages r < t csc , to for r > 2t csc . The fraction of light produced 
by stars inside clouds depends on the wavelength, the past star for- 
mation and chemical enrichment history, and the IMF. This fraction 
is much larger at shorter wavelengths (for example, the far-UV), for 
which the light is dominated by massive stars with short lifetimes. 
Since it is only for these stars that the extinction by clouds is signif- 
icant, we can approximate the recent SFR history when estimating 
this fraction. We assume the recent SFR to be constant in time for 
the case of a quiescent disk, while for an ongoing burst we use the 
actual exponential SFR history of that burst. The fraction of stars 
which are inside clouds are extincted both by the clouds and by the 
diffuse medium, while the remaining stars are extincted only by the 
diffuse medium. The extinction of stars inside clouds depends on 
the surface density of the clouds and on the gas metallicity, with 
the stars treated as being at the centre of each cloud. We then cal- 
culate the total dust attenuation by combining the attenuation fac- 
tors for the diffuse d ust and molecular clouds. For more details see 
lLacev et"al] d2010al) . Important features of our method are that the 
dust extinction varies self-consistently with other galaxy properties 
such as size, gas mass, and metallicity, and that young stars suffer 
more extinction than old stars, even at the same wavelength. 

Our new extinction model thus depends on the parameters 
/cloud, Mdoud, ^doud and t csc . In fact, the dependence on M c i oud 
and r c ioud is only through the combination M c i ou d/rci ou d which 
determines the surface density of a single cloud. Our default val- 
ues for thes e parameters are the same as we used in our GRASIL 
modelling in lBaugh et al.l u2005h . namely / c ioud = 0.25, M c ioud = 



10 Mq, r c i ou d = 16pc and t csc = lMyr. The value t csc was ad- 
justed in that paper to match the far-UV LF of LBGs at z ~ 3, 
w hile the values of the other three dust parameters were chosen 
in iGranato et ai] (2000) by comparison with observations of lo- 
cal galaxies. With these parameters, the extinction optical depth at 
1500A for stars at the centre of a cloud is around 90 for solar metal- 
licity, so the far-UV light from such stars is almost completely ex- 
tincted. The net attenuation by clouds will therefore be insensitive 
to the exact values of M c i ouc i and r c ioud- We investigate the effects 
of variations in / c i ou d and t csc in ^3.2.3 1 



2.3 Photoionization feedback. 

The reionization of the intergalactic medium (IGM) by a photoion- 
izing background suppresses the amount of gas cooling in small 
halos in two ways: the resulting IGM pressure inhibits collapse of 
gas into halos, and the radiative cooling time for gas inside halos 
is increased by the p hotoionizing b ackground. These effects were 
modelled in detail bv lBenson et al.1 d2002l) . and found to be reason- 
ably well approximated by a simple model in which gas cooling is 
completely suppressed in halos with circular velocities V c < Vcrit 
at redshifts z < «reion- The original Baugh05 and Bower06 mod- 
els assumed 2 re ion = 6 and Vcrit = 60km s _1 or 50km s" 1 re- 
spectively. This value for the reionization redshift was motivated 
by measurements of Gunn-Peterson absorption in quasar spectra, 
which imply that reionization is es sentially complete by z ~ 6 
dFan et alj[2000l : iBecker et alj|200lh . The values of the threshold 
halo circul ar veloc i ty, Vcr it, were motivated by the numerical sim- 
ulations of lGnedinl feoOO). Since then, measurements of the Thom- 
son scattering optical depth from fluctuations in the cosmic mi- 
crowave background h ave converged to imp ly a higher reionization 
redshift, z ~ 10 - 12 dDunklev et~al1 l2009). We therefore adopt a 
standard value z r cion = 10 in the current sudy. More recent numeri- 
cal simulations of the effect of reionization on collapse and cooling 
of gas in low mass ha los imply weaker photoionization feedback 
effect s than found by iGnedinl jHoeft et alj [20061: lOkamoto et al.l 
2008). In line with these simulations, we adopt Vcrit = 30km s _1 
as our standard choice in the present work. For brevity, we will 
refer to the slightly modified Baugh05 and Bower06 models with 
Zrcion = 10 and Vcrit = 30km s _1 simply as Baugh05 and 
Bower06 in what follows. We have not changed any other model 
parameters from their original values. 

2.4 Halo merger trees 

GALFORM calculates the evolution of galaxies in halo merger 
trees which describe the assembly and merger histories of dark 
matter halos. The halo merger trees are either computed using a 
Monte Carlo met hod based on the extended Press-Schechter model 
dCole e t al. 2000), o r extra cted fr om N-body simulati ons of the dark 
matter dHellv et al] 120031) . The iBaugh et ail d2005l) results were 
based on Monte Carlo trees, while those of iBower et all ( 2006) 
were ba sed on N-body trees extracted fr o m the Millennium sim- 
ulation dSpringel et al]|2005l) . ICole et alj d2008h found some dif- 
ferences between the Monte Ca rlo and N-body trees, as a result 
of which Parkinson et al. (2008) developed a modified version of 
the ICole et alj 12000) Monte Carlo algorithm, which brought the 
Monte Carlo trees into much better agreement with the N-body 
trees from the Millennium simulation. In this paper, we run both the 
Baugh 05 and Bower06 model s on Monte Carlo trees generated us- 
ing the lParkinson et al.1 g2008l) method. This has the advantage that 
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we can resolve much smaller halos at high redshifts than is possible 
using the Millennium simulation, which only resolves halos larger 
than 2 x 10 1 h~ Mq. In contrast, for the Monte Carlo trees used 
here, we resolve progenitor halos larger than 8 x 10 7 h" 1 Mq for 
model galaxies output at z — 3, 1.5 x 10 7 7i _1 Mq at z — 6, and 
4x 1O 6 /i _1 M at z = 10. This higher mass resolution is important 
for predictions of LBGs at high er redshifts and lower luminosities. 
The change to IParkinson et al \ Monte C arlo t rees leaves the pre- 
dictions published in Baug hetalJ d2005h and lBower etafl (2006) 
unchanged. 



3 EVOLUTION OF LBG LUMINOSITY FUNCTION 

In this section, we show the predictions of the models for the evo- 
lution of the rest-frame far-UV luminosity function (LF), and com- 
pare these with observational data for LBGs. The far-UV LF is one 
of the most basic observable properties of LBGs, since they are 
detected via their rest-frame far-UV emission, and is also of funda- 
mental importance, since the rest-frame far-UV emission is closely 
related to the SFR. For the purpose of this paper, we will define 
model LBGs to include all galaxies with far-UV (1500A) lumi- 
nosities in the relevant range, regardless of whether they satisfy the 
observational colour-selection criteria used in specific surveys, as 
discussed futher below. 

3.1 Comparison of Baugh05 and Bower06 models 

We start by showing in Fig. Q] the evolution of the dust-extincted 
rest-frame 1500A LFs over a large range in redshift and luminosity, 
for both the Baugh05 and Bower06 models (top and bottom panels 
respectively). For clarity, we show the evolution separately for 3 < 
z < 7 (left panels) and 7 < z < 20 (right panels). The solid lines 
show the total LF, while the dashed lines show the contribution 
to the LF from ongoing bursts. There are a number of interesting 
features in these plots. Both models show relatively mild evolution 
in the far-UV LF in the redshift range 3 < z < 8, with a gradual 
decrease in comoving number density with increasing redshift at 
most luminosities, followed by a much steeper decline in number 
density at z > 8. This rapid evolution at high redshifts is driven 
by the build-up of the dark matter halos hosting the LBGs, and 
was also found ealier in GALFORM predictions for the LFs of Lya- 
emitting gala xies (LAEs), which ar e a population closely related 
to the LBGs dLe Delliou et alfeOOrj) . However, the LFs in the two 
models have distinctly different shapes. 

In the Baugh05 model, the LF has a steep, quasi-exponential, 
cutoff at high luminosities, with a flattening or bump at lower lumi- 
nosities, followed by an upturn and steeper power-law slope at the 
faintest luminosities. This behaviour is due to the major role played 
by merger-driven starbursts in the Baugh05 model. As can be seen 
from the dashed lines in the figure, starbursts dominate the bright 
end of the LBG LF, while quiescent galaxies (with star formation 
occuring in the disk) dominate the faint end. The LF of the star- 
bursts actually declines at low luminosities, which is what causes 
the bump seen in the LF. The steep faint end is due to the quiescent 
population. In the Bower06 model, on the other hand, the LFs show 
a more power-law behaviour over the ranges plotted here, with a 
gradual steepening from low to high luminosity. There is a sharper 
break at the highest luminosities, but this is at much higher lumi- 
nosities and lower number densities than in the Baugh05 model. 
These differences are due to various effects. In the Bower06 model, 
ongoing starbursts do not dominate the LBG LF at any luminosity 



once dust extinction is included. However, due to the much shorter 
burst timescales in the Bower06 model, there is an important con- 
tribution to the far-UV LF from bursts which have recently termi- 
nated. In such cases, star formation has ceased, and the gas and dust 
associated with the burst have already been consumed or dispersed, 
but many of the massive stars produced during the burst are still 
shining in the far-UV. Such recently terminated bursts in fact domi- 
nate the bright end of the far-UV LF in the Bower06 model over the 
whole range of redshifts plotted in Fig. Q] and are responsible for 
the larger number of objects seen at the highest luminosities com- 
pared to the Baugh05 model. The dust extinction in these objects 
is low, as discussed further below. The shallower slope seen at the 
faintest luminosities in the Bower06 model is due to the stronger 
SN feedback. Even the LF for ongoing starbursts has a flatter shape 
in the Bower06 model, because bursts are triggered mainly by disk 
instabilities rather than galaxy mergers. 

We next compare the predicted far-UV LFs from the models 
with observational estimates of the LFs derived from samples of 
LBGs at redshifts z = 3 — 10. The observational selection is typi- 
cally based on two colours, one of which straddles the Lyman break 
at the target redshift, and the other which measures the spectral 
slope longwards of the break, and is used to exclude contaminants, 
principally lower redshift galaxies and galactic stars, which have 
redder colours than expected for a star forming galaxy at the target 
redshift. The observational LFs we plot already include corrections 
by the original authors for the completeness as a function of red- 
shift and luminosity, as well as for absorption by the IGM (where 
appropriate), but do not include corrections for dust extinction. We 
therefore compare our model LFs, including dust extinction, di- 
rectly with the LFs inferred observationally. Each LBG survey uses 
its own set of filters, so in practice they measure the far-UV LFs 
at slighly different effective rest-frame wavelengths, in the range 
1350-1700A. We plot our model LFs at a fixed rest-frame wave- 
length of 1500A, so, in principle, we should correct the luminosi- 
ties for this small wavelength difference when comparing with the 
observational data. However, the spectra of LBGs are observed to 
be fairly flat in this wavelength range, so this correction is small, 
and we neglect it here. It is also possible that the observed lumi- 
nosity functions are effectively missing some galaxies due to the 
colour selection, even after the completeness corrections applied 
by the original authors (e.g. due to the galaxies being too red in 
the far-UV due to dust extinction). Different completeness correc- 
tions have been applied in different surveys, and this presumably 
accounts for some of the differences seen in the inferred LFs. In 
contrast, our model LFs include all galaxies at a given 1500A lu- 
minosity, regardless of their far-UV colour. We defer to a future pa- 
per a detailed study of the effect of applying different LBG colour 
selections directly to the models, and the effect this has on the com- 
pleteness of the galaxy samples at far-UV wavelengths. The obser- 
vational LFs are all plotted for the same cosmology (Q m ~ 0.3, 
SIa = 0.7) as used in the Baugh05 model (correcting to this cos- 
mology if needed). The cosmology used in the Bower06 model is 
slightly different, but converting the observational data to it would 
have little effect on the comparison, so we ignore this correction 
here. 

We compare the 1500A LFs from the Baugh05 and Bower06 
models with the observed far-UV LFs in Figs.|2]and[3] Fig.[2]shows 
redshifts z — 3, 4, 5, 6 and Fig.[3]shows z = 7, 8, 10. The Baugh05 
model is shown in blue, and the Bower06 model in magenta. In 
both cases, we plot the dust-extincted LF as a solid line, and the 
unextincted LF as a dashed line. The Baugh05 model fits the ob- 
servational data well over the whole redshift range z — 3 — 10, 
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M JB (1500A) - 51og(h) M,„(1500A) - 51og(h) 

Figure 1. Predicted evolution of the rest-frame far-UV (1500A) luminosity function of LBGs in the Baugh05 model (top panels) and the Bower06 model 
(bottom panels). In both cases, the left panels show the evolution for 3 < z < 7 and the right panels for 7 < z < 20, with different redshifts shown in 
different colours, as indicated by the key. The dashed lines show the contribution from ongoing starbursts. All of the LFs plotted here include dust extinction. 



when we include dust extinction. This is remarkable when we con- 
sider that w e have not changed any model parameters from those 
published in lBaugh et alj d2005l) . apart from z rc ion and Van (as de- 
scribed in i|2j, whose effect is anyway quite small. We used the ob- 
se rved far-UV LBG L F at z — 3 to constrain the model parameters 
in lBaugh et al.l J2005f). but we did not use observations of LBGs at 
any other redshift. In contrast, the Bower06 model is in much worse 
agreement with the observed LBG LFs at almost all of the redshifts 
plotted, predicting far more very luminous galaxies than observed. 
The Bower06 model in its original form can therefore be excluded 
based on this observational data. 

By comparing the solid and dashed lines, we can assess the 
effect of dust extinction on the predicted LFs. For the Baugh05 
model, the effect of dust extinction is very large over the lumi- 
nosity range covered by the observational data, causing the LF to 
be shifted faintwards by ~ 1.5 — 2.5 mag at the bright end. In 
contrast, the effects of dust extinction are generally smaller in the 
Bower06 model, especially at the highest luminosities. An impor- 
tant factor in this is that the gas metallicity of the LBGs is typ- 
ically much larger in the Baugh05 than Bower06 model, leading 



to higher dust-to-gas ratios. The higher metallicities, in turn, re- 
sult from the top-heavy IMF assumed for starbursts in the Baugh05 
model, which leads to a higher yield of heavy elements from Type II 
supernovae. However, another important factor is the much shorter 
timescale for bursts in the Bower06 model (~ 0.001 Gyr or less as 
against ~ 0.01 — 0.1 Gyr in the Baugh05 model). A consequence 
of this is that the bright end of the far-UV LF in the Bower06 
model is dominated by bursts which have terminated within the 
last ~ 0.01 — 0.1 Gyr. For these objects, most of the massive 
stars formed during the burst are still shining, but the dust which 
shrouded the burst while it was ongoing has gone. As a result, they 
have large far-UV luminosities but very little dust extinction. 

Using a different semi-analytical model, Eo Faro et al. I d2009l) 
predicted a significant excess of faint LBGs compared to observa- 
tional data at z ~ 4 — 5. We do not see such an excess in our 
own model, even though we compare to similar observational data. 
This is most likely due to the dif ferent treatm ents of supernova and 
photoionization feedback in the lLo Faro et al J model. The effects of 
varying the treatment of supernova and photoionization feedback in 
our own model are discussed below in ^3.2. 1 l and i j3.2.5l 
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Figure 2. Predicted rest-frame 1500A LFs for the Baugh()5 model (blue) and Bower06 model (magenta) compared to observational data from LBG surveys. 
Solid and dashed lines show LFs with a nd without dust extin ction, (a) z = 3. (b) z = 4. (c) z = 5 . (d) z = 6. The observational d ata are as follows 
(with rest-frame wavelength): z = 3 -lArnouts et alj 120051) (crosses, 1500A). ISawicki & Thompson! 120061) (em pty triangles, 1700A). iReddv & S teidel 

4 -ISteidel et al " |l99^) (emptv circles, 1700A). ISawicki & Thompson! llOOfh (empty triangles, 17 00A).lYoshida et alj 

5 - lYoshida et alj j20O6l) (emp ty squares, lSOOAl | lwata et alj i2007h 
6 ^ IShimasaku et alj 12005) (open squares, 



2009) (filled circles. 1700A); .. 

20061) (empty squares, 1500A\lBouwens et al.H2007l) (fill ed circles, 1600A); 



(stars, 1 600A).lBouwens et alj i2007l) (filled circles, 1 600A).lMcLure et alj JiOQ^) (empty triangles, 1500A); 
1400A). lBouwens et alj <2007l) (filled circles, 1350Al lMcLure et alJfcOO^ (empty triangles, 1500A). 



3.2 Effects of varying model parameters 



3.2. 1 Photoionization feedback 



We now consider the effects on the far-UV LFs of varying some of 
the key parameters in GALFORM. Since the default Baugh05 model 
has been shown to agree much better with the observed LBG LFs 
than the default Bower06 model, we consider here only parameter 
variations around the default Baugh05 model. Our purpose here 
is to understand the sensitivity of the LBG predictions to different 
parameters. We emphasize that most of the parameter variations we 
consider do not lead to models that we would consider acceptable 
overall, since we require our model to fit a much wider range of 
observational data than just LBGs, as discussed in Jj2] 



Photoionization feedback is described by the parameters z r cion and 
Krit- In Fig. [4] we compare the far-UV LFs for our default val- 
ues of these parameters (z rc i on = 10, Krit = 30km s , blue 
curves) with the values used in the original iBaugh etaL (2005) 
paper (z rc ion = 6, V c -At = 60km s _1 , green curves), and with 
the case 2r e ion = 10, Krit = 60km s _1 (red curves). We only 
show this comparison at z — 3 and z = 6; for z ^ 10 the 
three models are identical. At z = 6, the model with z TC i on = 6 
and T4 r it = 60km s _1 is slightly above the default model at the 
faint end of the LF. This is because photoionization feedback has 
only just turned on in the former model, and so has not had time 
to have any effect on galaxy luminosities. On the other hand, by 
z = 3, the model with z rc i on = 6 and V C rit = 60km s _1 is 
somewhat below the default model at the faint end. This is be- 
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Figure 4. Predicted rest-frame 1500A LFs for the Baugh05 model, showing the effect of varying the photoionization feedback parameters, at (a) 2 = 3 and 
(b) 2 = 6. The blue lines show the default model (z re j on = 10, V cr i t = 30km s —1 ), green lines show the case where z rc j on = 6, V cr it = 60km s —1 , 
and red lines show the case z rc i on = 10, V C rit = 60kms _1 . The solid and dashed lines show the LFs with and without dust extinction, as before, and the 
observational data are the same as plotted in Fig. [2] 



cause the effect of photoionization feedback is much stronger with 
Vcrit, — 60km s~ than Vcrit = 30km s~ (the halo mass aff- 
fected scales approximately as V C rit/(l + z) 3 ' 2 ). The model with 
Zreion = 10 and Vcrit = 60km s -1 shows a large suppression of 
the faint end of the LF at z — 6 relative to the default model, but 
by 2 = 3 the predicted LF is very similar to that for 2reioi> = 6 and 
V^t = 60km s . This shows that the predicted LF becomes in- 
sensitive to the value of z le ion at much later ep ochs. The reduction 
in Vcrit for the present default model relative to lBaugh et al.l fc005h 
also causes a slight steepening in the present-day galaxy luminosity 
function, but only for galaxies fainter than Mb — 5 log h ~ —16. 
For the Bo wer06 model, the cha nge in Vcrit and 2 re ion from the val- 
ues used in iBower et alj J2006I) has negligible effects on either the 
LBG LFs shown here or on the 2 = luminosity functions. This is 
because the stronger SN feedback assumed in the Bower06 model 
dominates over the effect of photoionization feedback in low-mass 
galaxies. 

3.2.2 IMF and starbursts 

Starbursts and a top-heavy IMF play crucial roles in the Baugh05 
model, so we next consider the effects of changes in these com- 
ponents in Fig. [5] We show comparisons of the far-UV LFs for 
2 = 3, 6, 10 to span the range of current observational data on 
LBGs. In these panels, the blue line shows the default model with 
starbursts triggered by all major galaxy mergers and some minor 
mergers, with a top-heavy IMF in all of these starbursts. W e em- 
phasize that, acording to the results of Bau gh et al. | d2005h . all of 
these ingredients are required in order that the model reproduce the 
number counts and redshift distributions of sub-mm galaxies. The 
magenta line shows the effect of turning off starbursts completely. 
In this case, the predicted LFs are far below the observations. This 
is because, in the Baugh05 model, we deliberately chose a star for- 
mation timescale in quiescent disks which is long compared to the 
Hubble time at high redshift, so that disks at high redshift are gas- 
rich, providing more fuel for star formation in starbursts. As a less 
extreme variation, we consider allowing bursts to be triggered by 



major galaxy mergers only, shown by the green lines. In this case, 
the changes in the LF are smaller but still significant (factors of 
a few at all plotted luminosities), implying that most LBGs in the 
model are bursts triggered by minor mergers. Finally, we show by 
the red lines the effects of assuming a Kennicutt IMF in bursts (the 
same as for quiescent star formation), rather than a top-heavy IMF. 
In this case, the effects on the LF are mainly at the bright end. 
Compared to the default model, there is a large change in the unex- 
tincted LFs, but a smaller change after dust extinction is included. 
This is because the top-heavy IMF in bursts has two effects on the 
far-UV LF which partly cancel each other: there are more massive 
stars, and so higher intrinsic UV luminosities, but there is also more 
dust, due to increased metal production, and so more dust extinc- 
tion. The relatively modest differences between the dust-extincted 
LBG LFs with and without the top-heavy IMF show that fitting 
these data does not by itself provide a strong argument for intro- 
ducing variations in the IMF, although the fit is clearly better with 
the top-heavy burst IMF. The slope of the top-heavy IMF assumed 
in the Baugh05 model, x = (£[2j>, is quite extreme. However, the 
exact value of the slope is not important; any IMF with a low-mass 
cutoff > 5Mq would produce similar results. 

3.2.3 Dust extinction 

As already discussed, dust extinction has a large effect on the pre- 
dicted LBG LF. In Fig. [6] we show the effects of varying the two 
key adjustable parameters in our dust model, fdoud, the fraction 
of the ISM in molecular clouds (which affects the attenuation by 
the diffuse dust), and t csc , the timescale for stars to escape from 
the clouds in which they are born (which affects the attenuation by 
molecular clouds). The blue curves show the LFs for the default 
values (/doud = 0.25 and t csc = lMyr). The green and cyan 
curves show the effect of increasing t csc to 3 Myr and 10 Myr re- 
spectively, while keeping / c i ou d fixed at 0.25, while the red and 
magenta curves show the effect of increasing / c i ou d to 0.75 and 1 
respectively, while keeping t csc fixed at 1 Myr. The dust-extincted 
LF is seen to be quite sensitive to increasing the escape time from 
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Figure 3. Predicted rest-frame 1500A LFs for the Baugh05 model (blue) Figure 5. Predicted rest-frame 1500A LFs for the Baugh()5 model, showing 

and Bower06 model (magenta) compared to observational data for LBGs. the effects of varying the parameters for bursts, (a) z = 3, (b) z = 6. 

Solid and dashed lines show LFs with and without dust extinction, (a) (c) z = 10. Blue lines - default model; red - Kennicutt IMF in bursts; 

2 = 7. (b) z = 8. (c) z = 10. The observati onal data ar e as follows green - bursts triggered by major mergers only; magenta - no bursts. The 

(with rest-frame wavelength): 2 = 7- I puchi et all (2009) (empty cir- observational data are the same as plotted in Figs. [2]and[3] 

cles and uppe r limits , 1500A). IOesch et alj feOlObh (filled circles, 1600A) 
iMcLure et all feOld) (empty triangles, 1500 A); 2 = 8 - [ Bouwens et al J 
bOlCh (filled circles and upp er limits, 1700A) lMcLure et al.N20ld) (empty 

triangles, 1500A); 2 = 10 - iBouwens et al l (2009bJ) (filled circles and up- g, 00 oo RAS, MNRAS 000, 000-000 

per limits, 1600A). 
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Figure 6. Predicted rest-frame 1500A LFs for the BaughOS model, showing 
the effects of varying parameters controlling dust extinction in the model, 
(a) z = 3. (b) z = 6. (c) z = 10. The dashed blue line in each panel 
shows the LF without dust extinction, while the solid lines show LFs with 
dust extinction parameters as follows: blue - default model (t CBC = lMyr, 
/cloud = 0.25); green - t csc = 3Myr; cyan - t csc = lOMyr; red - 
/cloud = 0.75; magenta - / c i oud = 1. 



its default value, reflecting the fact that a larger fraction of the UV 
light is absorbed within the clouds when t csc is increased. This 
effect is especially strong in the Baugh05 model due to the top- 
heavy IMF assumed for bursts, which means that the UV emission 
is dominated by even higher mass stars than for a normal IMF. For 
the x — IMF, 90% of the 1500A light (integrated over time) is 
emitted by stars with lifetimes < lOMyr, as compared to lifetimes 
< lOOMyr for the Kennicutt IMF. The LF is also seen to be in- 
sensitive to modest increases in /cioud from its default value, but to 
become more sensitive as /cioud approaches 1. This is because the 
optical depth of the diffuse dust is proportional to 1 — /cioud- 

It can be seen that the effects of i OS c and / c i ou d on the far-UV 
LF are to some extent degenerate, with a decrease in t eac ha ving 
a similar effect to an increase in / c i oud . In lBaugh et al.l d2005t) . we 
chose to try to match the far-UV LFs of L BGs by reducing t e3C 
in bursts from the value 10 Myr assumed in lGranato et al.l ( 2000) 
to 1 Myr, while keeping the same value of /cioud = 0.25 as in the 
latter paper. We could probably obtain a similarly good fit to the ob- 
served LBG LFs with a larger value of t csc combined with a larger 
value of / c ioud, but we do not pursue this here. In future, it might 
be possible to distinguish the effects of t csc and /cioud by compar- 
ing with observed SEDs of LBGs over a broader wavelength range, 
since the two components of the dust extinction have different de- 
pendences on wavelength. 

Individual galaxies in the Baugh05 model selected by their ex- 
tincted 1500A luminosities are also predicted to have large far-UV 
extinctions, typically ~ 1 — 2 mag over the range of luminosity and 
redshift plotted in Figs.|2]and|3] though with a large scatter. There 
is a trend for the typical extinction to decrease somewhat with in- 
creasing redshift. Note that selecting galaxies by their extincted lu- 
minosities automatically biases against including the galaxies with 
the highest extinctions, which can cause the average extinction for 
such samples to differ from that of the underlying population as a 
whole. The values of dust extinction predicted by our model are 



similar to the observational estimates for LBGs by Bouwens et al 
(2009a) at z ~ 3, but somewhat larger at z ~ 6. iBouwens et al 



estimate extinctions from observed UV continuum slopes, assum- 
in g the same relation b etween UV slope and extinction as found 
bv lMeurer et alj (1999) for local UV-selected starbursts. However, 
the UV continuum slope is also sensitive to the age, metallicity and 
IMF of the stellar population, which could be different in high red- 
shift galaxies compared to local samples. Fu rthermore, eve n many 
local gal axies show large deviat ions from the lMeurer et alj relation 
(e.g|Bell 2002; B uat et al.1120101) . For these reasons, we defer a de- 
tailed comparison with observational estimates of dust extinction 
in LBGs to a future paper in which we directly compare predicted 
and observed far-UV colours in LBG samples. 

3.2.4 Burst timescale and duration 

Since starbursts play a very important role in the Baugh05 model, 
we investigate the effects on the LBG LFs of changing the burst 
timescale and duration. In Fig.|7]we show the effect of varying the 
parameters /d yn and r*buist,min which control the SFR timescale 
in bursts. The blue lines show the default Baugh05 model, which 
assumes / dyn = 50 and r»burst,min = 0.2 Gyr. The green and 
cyan lines show the effect of reducing r*burst,min to 0.02 Gyr and 
0.005 Gyr respectively, while keeping /d yn = 50. This is seen 
to have only a small effect on the dust-extincted LF at z — 3 
and z — 6, but to have a large effect at z = 10, where reducing 
T"*burst,min results in a larger number of bright LBGs. The larger 
sensitivity to r»bufst,mta at higher redshifts is because the spheroid 
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Figure 7. Predicted rest-frame 1500A LFs for the Baugh05 model, showing 
the effects of varying the star formation timescale in bursts, (a) z = 3. (b) 
z = 6. (c) z = 10. Blue lines - default model (T*turst tnln = 0.2Gyr, 
/dyn = 50); green - T„ burstimin = 0.02Gyr; cyan - T, burstimin = 
0.005 Gyr; red - r, burstjmin = 0.005Gyr, / dyn = 2. 
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Figure 8. Predicted rest-frame 1500A LF at z = 3 for the Baugh05 model, 
showing the effects of varying the burst duration as a multiple of its e- 
folding time. Blue lines - default (n T = 5); green - n T = 1; red - n T = 5. 



dynamical time Td yn (which enters in the burst timescale as shown 
in eq.lQ}) is typically shorter at high-z. The red line in Fig. [7] 
shows the effect of reducing both /d yn and TVburst.min to the values 
/dyn = 2 and T„burst,min = 0.005 Gyr assumed in the Bower06 
model. In this case, the effect on the dust-extincted far-UV LF is 
dramatic at all redshifts - there are many more very luminous LBGs 
and slightly fewer faint LBGs compared to the default Baugh05 
model. The effects of dust extinction on the LF are also much 
smaller, which is the main reason for the larger number of bright 
galaxies in the extincted LF. The reason for the smaller dust ex- 
tinction when the burst timescale is greatly reduced is that the stars 
that emit the 1500A light have typical lifetimes ~ 0.01 - 0.1 Gyr, 
longer than the burst e-folding time in this case. Therefore the dust 
associated with the burst has mostly been consumed or ejected 
while these stars are still shining. This is the same effect as dis- 
cussed earlier in connection with the Bower06 model. The differ- 
ence in burst timescales, rather than the top-heavy IMF, seems to be 
the main reason why the Baugh05 model is much more successful 
than the Bower06 model in explaining the observed LBG LFs. 

Fig.[8]shows the effect of changing the duration of bursts while 
keeping the star formation timescales and e-folding times fixed. 
The blue lines show our default model, in which the burst duration 
is n T = 3 e-folding times, while the green and red curves show the 
effect of decreasing this to n T = 1 or increasing it to n T — 5. The 
effects are seen to be quite small, showing that the predictions for 
the far-UV LF are insensitive to details of the burst duration, but 
are more sensitive to the star formation timescale. 



3.2.5 Supernova feedback 

Fig. [9] shows the effect of varying the supernova feedback param- 
eters Vkot and Qhot- We only show results for 2 = 3, since the 
changes for z = 6 and 10 are similar or smaller for the luminosity 
range plotted here. As before, the blue lines show the default model 
with Vhot = 300km s -1 and « no t = 2. The green and cyan lines 
respectively show the effect of reducing or increasing l/ n ot by a fac- 
tor 2, which causes the gas ejection rate for a given star formation 
rate to decrease or increase by a factor 4. We see that the LF in- 
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Figure 9. Predicted rest-frame 1500A LF at z = 3 for the Baugh05 model, 
showing the effects of varying the parameters for supernova feedback. Blue 



lines - default; green - V^, 



150km s 1 ; cyan - Vh t = 600km s 



red - a h ot = 1; magenta - a hot = 3. 



creases or decreases at the faint end as Vhot decreases or increases. 
The effect is modest for the dust-extincted LF, but larger for the 
unextincted LF, showing that the effects of changing feedback on 
the stellar luminosities are to some extent compensated by changes 
in the dust extinction. The red and magenta lines show the effect 
of decreasing ahot to 1 or increasing it to 3. A larger ahot means 
that the feedback increases more strongly as the galaxy circular ve- 
locity decreases. Again, the galaxy LF increases or decreases at the 
faint end in the sense expected, though by a modest amount. We 
conclude that the predicted LFs in the Baugh05 model are rela- 
tively insensitive to the supernova feedback adopted, in the range 
covered by current observational data. 



4 SIZES AND OTHER PHYSICAL PROPERTIES OF 
LBGS 

We now present predictions for other physical properties of LBGs 
as functions of their rest-frame UV luminosities. For galaxy sizes, 
we show some results for both the Baugh05 and Bower06 mod- 
els, but for the other properties, we only show the Baugh05 model, 
since this model is in much better agreement with observations of 
both LBG luminosity functions and sizes. We make a detailed com- 
parison here with observational data on the sizes, but for reasons 
of space we make only brief comparisons with observational con- 
straints on the other properties, deferring more detailed compar- 
isons to future papers. 



4.1 Sizes 

Following their luminosity functions, the half-light radii of LBGs 
in the rest-frame far-UV are their most directly observable physical 
property, if HST imaging is available. GALFORM predicts sizes for 
the disk and bul ge compon ents of each galaxy using the method 
detailed in lCole et al.1 {2000). In summary, disk sizes are calculated 
assuming conservation of angular momentum of the gas that cools 
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Figure 10. (a) Predicted rest-frame UV (1500A) half-light radii of LBGs as 
a function of rest-frame 1500A absolute magnitude for the Baugh05 model. 
Both radii and absolute magnitudes include dust extinction. The lines show 
median values, and the error bars show the 10-90% range. Different red- 
shifts are shown in different colours, as indicated by the key, and have been 
slightly horizontally offset from each other for clarity, (b) and (c): Com- 
parison of predicted rest-frame UV half-light radii for the Baugh0 5 model 
(blue) a nd Bower06 model ( magenta) with observa tional from Oesch et al. 
1 20 10a) (filled circles) and Bouwens et al ] fc004ah (crosses). Two lumi- 
nosity ranges are shown, 0.12 < L(UV) / L» < 0.3 in (b) and 0.3 < 
L(UV)/L, < 1 in (c), equivalent to -18.9 < A/ A b(1500A) - 
51ogh < -17.9 and -20.2 < A/ A b(1500A) - 51og/i < -18.9 re- 
spectively. In these two panels, the solid lines show the mean radius in that 
luminosity range, and the dashed lines show the medians and 10-90% per- 
centiles (Baugh05 model only). The observational data are the mean values, 
with the errorbars showing the error on the mean. 
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out of halos, while bulge sizes are calculated assuming energy con- 
servation when bulges form by mergers or by disk instabilities. The 
method allows for the gravity of the disk, bulge and dark halo, and 
includes the effects of halo contraction due to the gravity of the 
disk and bulge. Half-light radii in different bands are then calcu- 
lated by combining the predicted half-mass radii and luminosities 
for the disk and bulge components, assuming an exponential profile 
for the disk and an r 1//4 -law profile for the bulge. 

In the top panel of Fig. [10] we show the predicted median 
half-light radii (at rest-frame 1500A) as a function of rest-frame 
1500A luminosity for galaxies at different redshifts in the range 
z — 3 — 15, for the Baugh05 model. The stellar half-mass radii 
plotted in the same way are very similar to the half-light radii, so we 
do not show them here. (In principle they can differ if the disk and 
bulge components have different ages, metallicities or IMFs.) The 
galaxy sizes range from ~ kpc for the brightest LBGs to ~ lOOpc 
for the faintest and highest redshift LBGs. Both the Baugh05 and 
Bower06 models predict that sizes decrease with increasing red- 
shift at a given far-UV luminosity, although by factors that depend 
on the model and on the galaxy luminosity. This size evolution re- 
flects that for the host dark matter halos, for which the size scales as 
M^f o (l+z)~ 3 / 2 at the redshifts shown here, but it differs in detail. 
The ratio of galaxy to halo size i s not fixed, but d epends on the gas 
cooling and merger history (see lCole e t al. 2000 for more details). 
The far-UV luminosity is also only indirectly related to the halo 
mass, depending instead on the recent star formation rate, which 
may correlate only weakly with halo mass, especially if bursts are 
important, as in the Baugh05 model. 

The dependence of size on luminosity is distinctly different in 
the two models. In the Baugh05 model, it approaches r oc L 1 ^ 3 at 
the higher redshifts plotted, but is shallower at the lower redshifts. 
On the other hand, the Bower06 model shows a flat or even de- 
clining dependence of size on luminosity. These dependences are 
similar to what the respective models predict for the dependence of 
disk size on luminosity in the r-band at z = 0, where the Baugh05 
model was found to agree much better with the disk si ze-luminosity 
relation measured in the SDSS dGonzalez etai]|2009l) . 

In the lower two panels of Fig. [Toj we compare the size pre- 
dictions from both mod els with observational da ta on LBGs from 
iBouwens etal] ( 12004 al) and lOesch et all d2010ah . based on HST 
imaging in the rest-frame far-UV, and covering redshifts 2 < z < 
7. The observational data are given as mean half-light radii at dif- 
ferent redshifts for two different far-UV luminosity ranges, speci- 
fied as multiples of th e observed characteristic far-UV luminosity 
L* at z = 3 (taken bv lOesch et all to be Mab(1600A) = -21.0 
for h = 0.7). We present the model predictions in the same way. 
The two panels show the two luminosity ranges; the solid blue and 
magenta lines show the mean sizes predicted by the Baugh05 and 
Bower06 models respectively. The blue dashed lines show the me- 
dian sizes in the Baugh05 model, with the error bars indicating the 
10-90% ranges. For the higher luminosity range (0.3 < LjL* < 
1), both models are in reasonable agreement with the observations, 
but in the lower luminosity range (0.12 < L/L* < 0.3), only the 
Baugh05 model matches the observed sizes. This shortcoming of 
the Bower06 model results from the flat size-luminosity relation it 
predicts. The prediction by the Baugh05 model of LBG sizes in 
agreement with observations is a significant success of the model, 
since the only observed size information originally used in fixing 
the model parameters was disk sizes at 2 = 0. 

We now consider a range of other physical properties for the 
Baugh05 model only, shown as functions of dust-extincted rest- 
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Figure 12. Relation between Spitzer 3.6fim and optical (i or 2) apparent 
magnitude, for LBGs at 2 = 4 (red), 2 = 5 (green) and 2 = 6 (blue). 
The lines show the predicted median values in the Baugh05 model and the 
error bars show the 10-90% range. The filled circ les show the ob served 
median values (binned in optical magnitude) from I Stark et al ] 120091) . for 
the same redshifts as the model. For clarity, small horizontal offsets have 
been applied to lines for different redshifts, and the same offsets have been 
applied to the observational datapoints. 



frame far-UV luminosity for redshifts 3 < z < 15 in Figsll Hand 



us 



4.2 Stellar masses 

Median stellar masses are shown in the top left panel of Fig. QT| 
They cover a very wide range, from ~ 10 10 /i -1 iW© at the high- 
est luminosities and lowest redshifts, down to ~ 10 5 M© or 
less at low luminosities and high redshifts. The median stellar mass 
generally increases with increasing far-UV luminosity, but the de- 
pendence is flatter at lower luminosities and lower redshifts. At a 
given luminosity, the mass decreases with increasing redshift, but 
this dependence is much stronger at lower luminosities. 

In general, the stars in a galaxy are a mixture of populations 
formed quiescently and in bursts, with different IMFs. As we dis- 
cuss next, a large fraction of LBGs are predicted to be starbursts. 
In this case, the far-UV luminosity is generally dominated by the 
burst component, but the stellar mass need not be. For model LBGs 
at z — 3, ~ 50 — 80% of the stellar mass was formed quies- 
cently in disks, and even for those LBGs which are ongoing bursts, 
only 10-30% of the stellar mass has formed in the current burst. 
At higher redshifts and lower luminosities, the fraction of the stel- 
lar mass formed in bursts increases for galaxies selected by their 
far-UV luminosities. For example, at z = 6, the fraction of stel- 
lar mass formed quiescently increases with luminosity from 20% 
to 80% over the range plotted here, while for the subset which are 
ongoing bursts, the fraction formed in the current burst decreases 
from ~ 70% to ~ 10% over the same range. At z = 10, the stellar 
mass in LBGs is even more dominated by burst populations, with 
only ~ 10 — 30% formed quiescently, and ~ 50 — 90% formed in 
a current burst. 

The stellar masses of observed LBGs have been estimated in 
a number of studies, starting with z — 3 dSawicki & Yed[T998l : 
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Figure 11. Predicted physical properties of LBGs as functions of dust-extincted rest-frame 1500A absolute magnitude in the Baugh()5 model. Different 
redshifts are shown by different colours, as indicated by the key. Unless stated otherwise, all properties are plotted as medians, with error bars showing the 
10-90% range. The different panels are as follows: (a) stellar mass; (b) fraction of LBGs currently undergoing a burst; (c) dark matter halo mass; (d) mean 
bias (solid and dashed lines respectively show mean bias in a luminosity bin and for galaxies brighter than a given luminosity); (e) circular velocity of stars; 
(f) dust-extincted bulge-to-total luminosity ratio at 1500A (with bulge-to-total stellar mass ratio shown by dotted lines). 
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IPapovich et alj 2001: IShaplev et alj|2 001). and later extending to 
z = 4 - 6 dVerma et al.ll2007l:lEvles et alj|200l IStark et alj|2007l ; 
lYabe et al1 l2009: Star k et al.ll2009h . These various studies selected 
LBGs in different ranges of luminosity, using different colour se- 
lections, and obtained results which appear in conflict in some 
cases. However, all of them estimated stellar masses photometri- 
cally, by fitting models for galaxy SEDs (with an assumed IMF) 
to flux measurements at different wavelengths ranging, in the rest 
frame, from the far-UV to the optical. These SED models depend 
on a significant number of parameters in addition to the stellar 
mass, including dust extinction, age, star formation history and 
metallicity (and possibly redshift). In principle, all of these parame- 
ters should be determined simultaneously from the SED-fitting. In 
practice, this leads to degeneracies between different parameters, 
and sig nificant uncertainties i n the stellar masses, by factors up t o 
10 (e.g. lPapovich et alfeOOlUVerma et alJl2007l ; lYabe et alj|2009h . 
These studies all used a Salpeter or other similar Solar neighbour- 
hood IMF. If one allows larger variations in the IMF (as assumed 
in our model), then the uncertainties in the estimated stellar masses 
become even larger. Since the top-heavy IMF plays an important 
role in our model, it is not meaningful to compare stellar masses 
from the model directly with observational estimates which assume 
a Salpeter or similar IMF. Therefore we do not make a detailed 
comparison with values of stellar mass es timated from observa- 
tion al data, but simply no te a few numbers. IPapovich et all feOOlh 
and IShaplev et alj ( 1200 lh find typical stellar masses ~ 10 10 M© 
for LBGs a t z = 3 with M AB (1500A) - 5 log h ~ -20. 
IVerma et alj d2007l) find masses ~ 2 x 10 9 Mq for similar lumi- 
nosities at z = 5. For LBGs at z = 4 - 6. IStark et all 120091) 
find stellar masses ~ 10 9 - 1O 1O M over the luminosity range 
-19 < M A b(1500A) - 5 log h < -21, with only a weak depen- 
dence on redshift at a given luminosity. Our model predicts values 
of the stellar mass which are ~ 30 times smaller at a given lumi- 
nosity and redshift. 

However, the observational estimates of stellar mass are driven 
mainly by the fluxes measured at rest-frame optical wavelengths, so 
given the effects of the IMF, it is more meaningful to compare di- 
rectly with these flux measurements. As an example, in Fig.ll2lwe 
compare our model predicti ons with the Spitzer 3.6/im fluxes mea- 
sured bv lStark et al.1 ((2009) for LBGs at z = 4 — 6 as a function 
of the optical flux. We have chosen the lStark et alj sample because 
it is large, homogeneously selected, and covers a wider range of 
redshift and luminosity than other studies. At the redshifts of these 
galaxies, the 3.6/im band corresponds to a rest-frame wavelength 
of 5000 — 7000A, while the i- or z-band optical magnitude cor- 
responds to 1200 — 1600 A in the rest frame. In the Figure, the 
lines show the medians predicted by the model, while the filled 
circles show the median values of 3.6/im magnitude measured by 
IStark et all The predicted relation is seen to lie 0.8-1.6 mag (i.e. a 
factor 2-4) fainter in the 3.6 /im magnitude than the observations. 
The discrepancy is thus ~ 10 times smaller than the apparent dis- 
crepancy in stellar masses, confirming that the latter is mostly due 
to the difference between the IMF in the model and the IMF as- 
sumed in the observational estimates of stellar mass. Interestingly, 
the model predicts a very similar trend of rest-frame optical vs. far- 
UV luminosity as seen in the observational data, with almost no 
dependence on redshift in the range z — 4 — 6. The weak depen- 
dence on redshift in the model presumably results from the LBGs 
mostly being starbursts. We will make a more detailed comparison 
of the predicted rest-frame far-UV to optical SEDs of LBGs with 
observations in a future paper. 



4.3 Burst fraction 

The top right panel of Fig. II llshows the fraction of galaxies which 
are currently undergoing a burst. The burst population dominates at 
higher luminosities and higher redshifts, with a transition to more 
quiescent systems at lower luminosities. This is important for un- 
derstanding many of the other properties. In the Baugh05 model, 
bursts are triggered only by galaxy mergers. As discussed in i]3.2.2| 
the bursts responsible for LBGs in the model are typically triggered 
by minor, rather than major, galaxy mergers. The burst e-folding 
timescale varies over the range ~ 0.01 — 0.2Gyr for model LBGs 
in the luminosity range shown here, for z — 3 — 10. It typically 
increases with luminosity but decreases with increasing redshift. 
Note that the burst e-folding time can be significantly shorter than 
the burst SFR timescale due to the effect of supernova feedback. 
Bursts are assumed to last for 3 e-folding times in our fiducial 
model. However, the brightest LBGs are on average seen only a 
small fraction of an e-folding time after the burst began - this is 
presumably a selection effect due to a burst being brightest in its 
early stages when the SFR is highest. 



4.4 Halo masses and clustering bias 

Dark halo masses are shown in the middle left panel of Fig. QT] 
The trends with luminosity and redshift are similar to those already 
discussed for the stellar mass. They also cover a wide range, from ~ 
1O 12 /i _1 M0 at the highest luminosities and lower redshifts, down 
to ~ 10 9 7i _1 Mq or less at low luminosities and high redshifts. 

The masses of the dark halos hosting LBGs are not directly 
observable, but they can be constrained from clustering measure- 
ments. In the middle right panel of Fig.QT| we show the predicted 
mean large-scale clustering bias of LBGs. We ca l culate d this using 
the analytical halo bias formula of Isheth et alj |2001), averaging 
the bias over the distrib ution of host halo masses in each bin or 
range of luminosity (see lBaugh et alj fl998 for details). This bias 
applies on scales larger than roughly the sum of the virial radii 
of the two host halos. The solid lines show the mean bias within 
a luminosity bin, while the dashed lines show the mean bias for 
galaxies brighter than a given luminosity. The bias increases with 
increasing redshift, and also generally with increasing luminosity, 
although the latter dependence is weak at lower redshifts. A similar 
behaviour for bia s of LAEs in the same model was earlier found by 
lOrsi et al We predict a bias b ~ 2 at z = 3, nearly inde- 

pendent of luminosity, increasing to b ~ 4 — 7 at z — 10 and to 
b ~ 8 at z = 15. This increase of the bias with redshift reflects 
the fact that selecting at a fixed luminosity picks out host halos fur- 
ther and further out on the tail of the halo mass function at larger 
and larger redshifts, which are more and more strongly clustered 
relative to the dark matter as a whole. 

Obs ervational measurement s of c lustering of LBGs incl ude 
those of lAdelberger et al. (1998) and iGiavalisco etal.1 dl998l) at 
z = 3 and lOuchi et alj j2004h" at z — 4 — 5. IOuchi et alj combine 
their own and previous clustering measurements with theoretical 
predictions of the dark matter clustering to estimate b ~ 2.7 ± 0.4 
at z » 3, b ss 3.5 ± 0.7 at z w 4 and b « 4.6 ± 1.1 at z « 5 
for LBGs with M AB (1500 A) - 5 log ft < -20. Our model predic- 
tions are b ~ 2.2, 2.8, 3.5 respectively for the same redshifts and 
luminosities, which are consistent with these estimates within their 
errors. We will investigate the clustering in more detail in a future 
paper. 
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4.5 Circular velocity 

The circular velocity at the stellar half-mass radius is plotted in 
the lower left panel of Fig. QT| While the trends with luminosity 
and redshift are qualitatively similar to those for stellar and halo 
mass, the dependence on redshift is weaker and the total range 
of values is much smaller, from ~ 40km s -1 at the lowest lumi- 
nosities and highest redshifts plotted, up to ~ 250km s -1 at the 
highest luminosities and lowest redshifts. At higher luminosities 
(Mab(1500A) — 5 log h < —19), there is almost no dependence 
of circular velocity on redshift. 

The velocity widths of LBGs at z ~ 3 have been measured 
from their rest-frame optical emission lines bv lPettini et al.l < fl~998L 
1200 lh . These lines (unlike emission and absorption lines in the rest- 
frame UV) reflect the kinematics of the star-forming gas in the 
galaxie s, which should be related to the rotation velocity of the 
galaxy. IPettini et al] teOOll) measure a median 1-D velocity disper- 
sion a » 70kms _1 in LBGs with Af AB (1500A)-5 log ft ~ -21. 
Allowing for inclination and other effects in galaxy disks , this cor- 
respo nds to a typical circular velocity V c ~ 120km s _1 dRix et all 
1 19971) . The line widths measured from CO in t wo gravitationalry- 
lensed LBGs imply similar values for a and V c dBaker et"al1l 20041 ; 
Coppin et al. 2007). Our model predicts a median circular veloc- 
ity ~ 200km s" 1 at the same luminosity and redshift, somewhat 
larger than the observational estimates. 

4.6 B/T and morphology 

The bottom right panel of Fig. [TT] shows the bulge-to-total ratio 
B/T, either measured in dust-extincted rest-frame 1500A luminos- 
ity (solid lines) or in stellar mass (dotted lines). We see that LBGs 
are very bulge-dominated in their rest-frame UV light at higher lu- 
minosities, but with a transition to disk domination at low luminosi- 
ties. This reflects the dominance of bursts at the higher luminosities, 
since our model assumes that star formation in bursts all occurs in 
the bulge component. The dominance of the bulge is generally less 
extreme in terms of stellar mass. At the highest luminosities, LBGs 
are predicted to be disk-dominated in stellar mass even though they 
are bulge-dominated in rest-frame UV light. This reflects the trig- 
gering of many starbursts by minor galaxy mergers, which leave 
the stellar disk of the larger galaxy intact, but sweep all of the cold 
gas into the bulge, where it forms stars in a burst. 

H ST imaging has revealed a wide range of morphologies for 
LBGs. Ravindranath et al. (2006) fit Sersic profiles to the rest- 
frame far-UV images of a large sample of LBGs at z = 3 — 5 with 
Mab(1500A) - 5 log h < -19.5. Based on this, they classified 
about 30% as bulge-dominated, 40% as e xponential, and 3 0% as 
multiple cores, suggesting galaxy mergers. lLotz et alj J2OO6I) anal- 
ysed a smaller sample of LBGs at z ~ 4 (also in the rest-frame far- 
UV) using completely different techniques, but arrived at similar 
conclusions, finding ~ 30% undisturbed bulgelike morphologies, 
~ 10 — 25% major mergers, and ~ 50% exponential disks or minor 
mergers. Our model predicts that most LBGs at these luminosities 
should be minor and major mergers, with the far-UV light domi- 
nated by a bulgelike burst component, which seems qualitatively 
consistent with the observational results. 

4.7 Star formation rates and specific SFRs 

We show the SFRs of LBGs in the top left panel of Fig. [13] The 
solid lines show SFR plotted against dust-extincted far-UV lu- 
minosity, and the dotted lines against the unextincted luminosity. 



There is an almost linear relation between the SFR and the far-UV 
luminosity, with the lines for different redshifts lying almost on top 
of each other. A constant linear relation between rest-frame 1500A 
luminosity and SFR would be expected under the following con- 
ditions: (a) a single IMF dominates; (b) dust extinction is zero or 
constant; (c) SFR varies on timescales > 10 8 yr. In the Baugh05 
model, stars form with a top-heavy (x = 0) IMF in bursts and a 
Kennicutt IMF in quiescent disks. For a constant SFR and Solar 
metallicity, the ratio L U (1500A)/SFR is 3.4 times larger for the 
x — IMF than for the Kennicutt IMF (and 2.6 times larger for 0.2 
times Solar metallicity). The LBG LF changes from being domi- 
nated by bursts at high luminosities to being dominated by quies- 
cent galaxies at low luminosities, and this causes a corresponding 
shift in the unextincted SFR vs L„(1500A) relation, which can 
be seen for the lower redshifts in the plot. The effects of changes 
in metallicity with redshift on the unextincted relation seem to be 
small. Comparing the solid and dotted lines, dust extinction is seen 
to introduce more scatter into the SFR vs L„(1500A) relation, but 
the average effect of the extinction depends only modestly on lumi- 
nosity and redshift. Observationally, SFRs of LBGs are generally 
inferred directly from their rest-frame far-UV luminosities, with or 
without a correction for dust extinction, so such estimates do not 
provide an independent test of the relation plotted here. Observa- 
tional studies typically assume a Salpeter IMF over the mass range 
0.1 — IOOA/0 for converting luminosities to SFRs. This Salpeter 
IMF would require an SFR 4.6 times larger than our top-heavy burst 
IMF to produce the same unextincted far-UV luminosity (for a con- 
stant SFR and Solar metallicity). 

The top right panel of Fig. Q~3] shows the specific star forma- 
tion rate (SSFR), defined here as the ratio of the current SFR to the 
stellar mass. There is a strong increase of the SSFR with redshift, 
and a weaker trend with luminosity. LBGs at higher redshifts are 
thus forming stars at a much larger fractional rate than lower red- 
shift LBGs, by factors up to ~ 10 3 . The age, t(z), of the universe 
is, of course, much less at higher redshift (shrinking from 2. 1 Gyr 
at z = 3 to 0.26 Gyr at z = 15), so it might be more physically 
meaningful to plot t(z) x SFR/M a t ar - This would still increase 
by a factor ~ 10 2 from z — 3 to z — 15. Finally, we should 
account for the fact that the rate of buildup of mass in long-lived 
stars is actually (1 — R) x SFR, where R is the fraction of the 
initial stellar mass returned to the ISM by mass loss from dying 
stars. R has the value 0.41 for the Kennicutt IMF but 0.91 for the 
top-heavy IMF. So the current rate of buildup of stellar mass com- 
pared to the past average rate is (l — R)t(z)SFR/M s tar- Allowing 
for the shift between burst and quiescent domination with chang- 
ing luminosity, this latter quantity still increases by a factor ~ 10 
over the range z — 3 to z — 15, with the high-z LBGs having 
(1 — R)t(z)SFR/M sta r ~ 10. Thus, by any reasonable measure, 
the highest redshift LBGs are predicted to be assembling stars ex- 
tremely rapidly compared to their past average rate. 

4.8 Cold gas masses and gas fractions 

The middle left panel of Fig. [T3] shows the cold gas masses of 
LBGs, where by "cold" gas we mean all of the gas which has con- 
densed into the galaxy, as distinct from the "hot" gas which remains 
distributed in the halo. Most of the cold gas will be in either atomic 
or molecular form. The brightest LBGs are predicted to have cold 
gas masses ~ 10 10 /i _1 Mq. At higher luminosities and higher red- 
shifts, there is a very nearly linear relation between the cold gas 
mass and the dust-extincted far-UV luminosity. This is a conse- 
quence of two effects: the nearly linear relation between SFR and 
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Figure 13. Further predicted physical properties of LBGs as functions of dust-extincted rest-frame 1500A absolute magnitude in the Baugh05 model. The 
lines show medians and the error bars indicate the 10-90% range. The different panels are as follows: (a) SFR - the solid and dotted lines show the relations 
with and without dust extinction in the UV luminosity; (b) specific SFR; (c) cold gas mass; (d) cold gas fraction; (e) gas metallicity; (f) stellar metallicity. 
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dust-extincted luminosity already noted above, and the fact that the 
more luminous LBGs are bursts, for which SFR = M ga s/r*, with 
most of the bursts having very similar star formation timescales, 
t*, according to our model. The linear relation between M gas 
and L„(1500A) breaks down at low luminosities where quiescent 
galaxies become important. These have longer SFR timescales than 
bursts, and so must have more gas to produce the same far-UV lu- 
minosity from young stars. 

The middle right panel of Fig. Q3] shows the gas fractions in 
LBGs, where we define this fraction as the ratio of cold gas mass to 
total cold gas + stellar mass in a galaxy. Although the gas fraction 
shows a very large scatter at lower luminosities, the median value 
is predicted to be very high ~ 90 — 99%. This results from two 
effects in the model: the SFR timescale in disks is large compared 
to the age of the universe at high- z, so that these disks are gas-rich, 
and most of the bursts are triggered by minor mergers, for which we 
require the gas fraction in the disk of the primary galaxy to exceed 
75% for a burst to be triggered. 

Molecular gas has been observed through its CO emission in 
two gravitationally lensed LBGs at z ~ 3, both of whi ch have rest- 
frame luminosities -Mab ( 1 500 A) - 5 log h ~ -20 (Bak er et ail 
2004!; ICoppin et alj |2007). The two galaxies have CO luminosities 
which differ by a factor 7. The conversion fro m CO lumino sity to 
molecular gas mass is significantly uncertain. Coppi rTet al.l use a 
conversion factor estimated for local starburst galaxies, and find 
molecular gas masses of 0.34 and 2.5 x 10 9 M Q for the two galax- 
ies. Using instead the conversion factor estimated for the Milky 
Way would increase these estimated gas masses by a factor ~ 6. 
Our model predicts a median cold gas mass (including both molec- 
ular and atomic gas) ~ 4 x 1O 9 M0 at the same luminosity and red- 
shift. This seems compatible with the current observations, given 
their uncertainties. 

4.9 Metallicities of gas and stars 

Finally, the bottom left and right panels of Fig. [T3] show the metal- 
licities of cold gas and stars respectively (the latter being a mass- 
weighted mean value over stars of all ages). We see that, at a 
given redshift and luminosity, the metallicity of the gas is gen- 
erally somewhat higher than that of the stars for LBGs that are 
dominated by bursts. This reflects the buildup of gas metallicity 
by self-enrichment in a current burst, while the stellar metallicity is 
an average over past activity. For both gas and stars, the metallici- 
ties typically depend only weakly on luminosity, and more strongly 
on redshift, again for the LBGs which are dominated by bursts. 
However, even at very high redshifts, the metallicities of LBGs are 
predicted to be non-negligible (e.g. Z gas ~ 0.003 at z — 15), due 
to self-enrichment by bursts which have large heavy element yields 
due to the top-heavy IMF. 

The gas metallic ities of z ~ 3 LBGs have been estimated 
bv lPettini et alj J2001) from rest-frame optical HII region emission 
lines, for a sample of 4 LBGs with A/ab(1500A) — 5 log ft ~ 
—21. They find significant uncertainties in the values for indi- 
vidual galaxies, but they are all constrained to be in the range 
~ 0.1 — IZq (where the solar metallicity Zq — 0.02), with apref- 
erence for somewhat subsolar values. Our model predicts roughly 
solar metallicities for the gas at this luminosity and redshift, com- 
patible with the upper end of the observed range. The stellar metal- 
licities of LBGs are only very weakly constrained by SED-fitting 
to observ ed broadband fluxes, due to degeneracies with other pa- 
rameters ( Papo vich et alj200ll) . but are compatible with our model 
predictions. 
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Figure 14. Predicted evolution of the LBG luminosity function at z = 7 — 
20 in the Baugh05 model, shown as a surface density of objects (number 
per solid angle per unit redshift) vs dust-extincted apparent magnitude at a 
fixed rest-frame wavelength of 1500A. 

5 PREDICTIONS FOR LBGS AT VERY HIGH 
REDSHIFTS 

In this section, we present additional predictions for LBGs at very 
high redshifts, z > 7, which are starting to be probed with HST, 
and which in future will be probed to higher redshifts and fainter 
fluxes by JWST, as well as by future 20-40m ELTs on the ground, 
including the E-ELT, TMT and GMT. NIRCam on JWST will do 
imaging over the wavelength range 0.6-5/im, and so can in princi- 
ple detect LBGs over the whole redshift range z = 7 — 20 and be- 
yond, if they are bright enough and numerous enough. Among the 
surveys planned with NIRCam are the Deep-Wide Survey (DWS), 
to reach tuab = 30 over an area of 100 arcmin 2 , and the Ultra- 
Deep Imaging Surve y (UDS), to reach ?tiab = 31 over 10 arcmin 2 
dGardner et alj2006h . 

We start by presenting in Fig. [14] predictions from the 
Baugh05 model for the number of z > 7 LBGs, this time shown as 
the surface density of objects (number per unit redshift per square 
arcmin) vs apparent magnitude. For simplicity, we present results 
in this section in terms of the apparent magnitude for a fixed rest- 
frame wavelength of 1500A, regardless of redshift. In reality, future 
LBG surveys will use a variety of filters, depending both on the 
telescope and on the redshift being targetted. However, in order to 
detect the Lyman-break feature reliably, such surveys will always 
include at least one filter probing wavelengths close to 1500A in 
the rest frame. Predictions for the number of LBGs in any future 
survey can therefore be approximately read off from Fig. [14] us- 
ing the apparent magnitude limit for the filter which is closest to 
rest-frame 1500A at the target redshift. We see from Fig. 1 141 that 
the planned DWS and UDS surveys on JWST should each detect a 
few LBGs at z = 15, and many more at lower redshifts. For con- 
venience in what follows, we note that for the cosmology assumed 
in the Baugh05 model, the relationship between apparent and ab- 
solute magnitudes at the same rest-frame wavelength is ulab ~ 
(Mab - 5 log h) = (46.16,46.37,46.71,46.97,47.29,47.69) at 
the redshifts z = (7, 8, 10, 12, 15, 20) plotted in Figs. [141 and [TBI 
An apparent magnitude uiab = 31 at z = 15 thus corresponds to 
an absolute magnitude Mab(1500A) — 5 log h = —16.3. 
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We next present in Fig. \T5\ a few directly observable prop- 
erties of high-jz LBGs as functions of rest-frame 1500A apparent 
magnitude. These particular properties are chosen because they are 
directly relevant to the detectability of these LBGs by different in- 
struments on the JWST and on future ELTs. The top panel of Fig.l 151 
shows the angular half-light radius at a rest-frame wavelength of 
1500A. This shows that at a given apparent magnitude, LBGs have 
slightly larger angular sizes at higher redshifts. This reversal of the 
trend seen for proper size vs absolute magnitude shown in Fig. 1 101 
is due to the combined effects of angular diameter and luminosity 
distances. LBGs at thab = 30 and z — 15 are predicted to have 
angular radii ~ 0.05 arcsec, close to the diffraction limit of JWST at 
the relevant wavelength. The middle panel shows the predicted cir- 
cular velocity at the stellar half-mass radius. This affects the widths 
of emission and absorption features in a galaxy spectrum, and so is 
relevant for spectroscopic studies. The dependence of circular ve- 
locity on apparent magnitude and redshift is again modest, varying 
from ~ 40km s _1 to ~ 150km s" 1 over the whole range plotted. 
At a fixed apparent magnitude, LBGs have larger circular velocities 
at high redshift. 

Finally, in the bottom panel we show the predicted Lya line 
flux. This depends on the fraction of Lya photons which escape 
from galaxies, which is very difficult to calculate from first prin- 
ciples due to the effects of resonant scattering of Lya by atomic 
hydrogen. We therefore follow the same approach as in our previ- 
ous work on Lya emitters using GALFORM, and assume a constant 
escape fraction, tuned to match observation s at lower redshifts. 
iLe Delliou etal] d2006h and lOrsi et ai] d2008h showed that the ob- 
served numbers and properties of Lya emitters at z = 3 — 6 could 
be reproduced well in the Baugh05 model by assuming a constant 
escape fraction of 2%, so we use the same value here. We see that 
LBGs at different redshifts all fall on the same linear relation be- 
tween Lya flux and rest-frame far-UV flux. LBGs with rrihs, = 30 
are predicted to have Lya fluxes ~ 10 _19 erg cm~ 2 s~ 1 . This is 
probably too faint to be detectable with JWST but should be within 
the reach of the E-ELT. 
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Figure 15. Predicted observable properties of LBGs at 2 7 in the 
Baugh05 model as functions of dust-extincted apparent magnitude at a fixed 
rest-frame wavelength of 1500A. The lines show median values and the er- 
ror bars show the 10-90% ranges, (a) Angular half-light radius at rest-frame 
1500A. (b) Circular velocity at stellar half-mass radius, (c) Lya flux (as- 
suming a Lya escape fraction / csc = 0.02). 



6 EVOLUTION OF SFR AND UV LUMINOSITY 
DENSITIES 

The final topic we consider is the cosmic star formation history and 
how this is traced by the far-UV luminosity density. We show in 
Fig.[T6]the comoving SFR density as a function of redshift for the 
Baugh05 and Bower06 models (top and bottom panels). The solid 
blue curves show the total SFR in all galaxies, while the solid green 
and red curves show the separate contributions from star formation 
in disks and in bursts. 

In the Baugh05 model, the total SFR density peaks at z ~ 3, 
with quiescent galaxies dominating at z < 2.5 and bursts dominat- 
ing above this. Similar results were shown in[Baugh et ai] d2005h 
and lLacev et al. I l l2010th . but there are small differences of detail 
because of the modified values for the photoionization feedback 
parameters, z re ion and V ^ rit, tha t we have used in this paper. As 
discussed in lLacev et al.1 feOlObl) . we also find it useful to show 
the SFR density in massive stars only (dashed lines), which we 
define as stars with masses m > 5Mq. These stars, which have 
lifetimes < 1 x 10 8 yr, dominate the production of UV radiation. 
The two IMFs in our model (assumed to cover the stellar mass 
range 0.15 < m < I2OA/0) have very different fractions of their 
initial stellar mass in high mass stars: f(m > 5Mq) — 0.24 
for the Kennicutt IMF assumed for quiescent star formation, and 
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Figure 16. Comoving SFR density vs redshift in the Baugh05 and Bower06 
models (top and bottom panels). In each panel, the blue line shows the total, 
while the green and red lines show the separate contributions of disks and 
bursts. The solid lines show total SFRs integrated over all stellar masses, 
while the dashed lines show the SFRs in high-mass (m > 5Mq) stars 
only. 

f(m > 5M0) = 0.96 for the top-heavy IMF assumed for bursts. 
The SFR density for massive stars evolves more strongly than that 
for all stars, increasing by a factor ~ 20 from 2 = to its peak at 
2 ~ 3, and then declining by a factor ~ 100 to 2 = 15. 

In the Bower06 model, the total SFR density peaks at almost 
the same redshift as in the Baugh05 model, though at a slightly 
lower value. However, unlike in the Baugh05 model, the SFR den- 
sity is dominated by quiescent star formation at all redshifts, and 
it also drops off more steeply with redshift beyond the peak. The 
SFR density for massive (m > 5Mq) stars in the Bower06 model 
is roughly a factor 3 lower than in the Baugh05 model for 2 < 5, 
but this difference increases at higher redshifts, reaching a factor 
~ 20 at 2 = 10. 

In Fig. [T7] we show the predicted evolution of the rest-frame 
1500A emissivity (i.e. luminosity density), for the Baugh05 model 
only. In the left panel, we show the luminosity densities with and 
without dust extinction (solid and dashed lines), and also show the 
separate contributions of quiescent and bursting galaxies to the to- 



tal. Dust extinction has a large effect on the far-UV emissivity, with 
the mean extinction (defined from the ratio of emissivities with and 
without dust) increasing from ~ 1 mag at 2 = to w 2 mag at 
2 = 3 — 15. The mean extinction in the bursts is larger than that in 
quiescent galaxies at most redshifts, but at 2 > 1, the bursts dom- 
inate the far-UV emissivity with or without dust extinction. As ex- 
pected, the unextincted far-UV emissivity approximately tracks the 
SFR density in high-mass stars shown in Fig.[l6](although changes 
in the metallicity and in the mix of the IMFs with redshift mean that 
the scaling is not exact). For comparison, in the Bower06 model, 
the unextincted far-UV emissivity is below that for the Baugh05 
model at all redshifts. However, the mean dust extinction is also 
lower, peaking at ~ 1 mag at 2 ~ 3, and falling to ~ 0.3 mag 
at 2 = 10. As a result, the extincted emissivity is similar to the 
Baugh05 model at 2 < 7, but falls below it at higher redshifts. 

In the right panel of Fig. [T7] we show the contributions 
to the dust-extincted 1500A emissivity from galaxies with rest- 
frame 1500A absolute magnitudes brighter than A4ab(1500A) — 
5 log ft = —14, —16 or —18, again for the Baugh05 model. An 
LBG survey reaching down to A/ab(1500A) —5 log ft < —18 (the 
typical limit for current surveys at 2 > 7) will detect the galaxies 
responsible for 50% of the far-UV emissivity at 2 = 5, but only 5% 
at z — 10. Detecting 80% of the far-UV emissivity at 2 = 10 re- 
quires detecting galaxies down to Mab( 1500 A) — 5 log h < —14, 
corresponding to an apparent magnitude ttiab ~ 33 at the same 
rest-frame wavelength. 

Finally, we consider the production of hydrogen-ionizing Ly- 
man continuum (Lyc) photons, which is of critical importance for 
reionizing the IGM. Fig. [T8] shows the evolution of the Lyc emis- 
sivity from galaxies predicted by the Baugh05 model. We plot the 
emissivity without applying any correction for absorption of Lyc 
photons by dust or gas in the galaxies, (i.e. assuming an escape frac- 
tion of 100%). The actual escape fraction is probably dominated 
by absorption on neutral hydrogen, but is currently very uncertain 
both theoretically and observationally (e .g.[Loeb & B arkanalboOll ; 
iBenson et al"]|2006l : Ishaplev et alj|2006h . The left panel shows the 
Lyc emissivities of quiescent and bursting galaxies as green and 
red lines, with the total in blue. The Lyc emissivity is even more 
dominated by bursts than is the 1500A emissivity (at least if ab- 
sorption by dust and gas are ignored). Bursts dominate at all red- 
shifts 2 > 0.2, and the ionizing emissivity of bursts is 7-50 
times larger than that of quiescent galaxies over the whole range 
2 = 3 — 15. The reason for this is that stars formed with the top- 
heavy x = IMF produce 1 1 times more Lyc photons per unit 
mass of stars formed than the Kennicutt IMF (for Solar metallic- 
ity). For comparison, the ratio is only 3.4 for production of 1500A 
photons. The right panel shows the contributions to the Lyc emis- 
sivity from galaxies brighter than various absolute magnitude limits 
in dust-extincted 1500A light. For example, galaxies brighter than 
M A b(1500A) - 5 log ft < -18 are predicted to emit 40% of the 
hydrogen-ionizing photons at 2 = 5, but only 3% at 2 = 10. Re- 
solving 80% of the ionizing emissivity at 2 = 10 requires detecting 
galaxies down to A/ab(1500A) — 5 log ft < —14 (corresponding 
to ttiab ~ 33), as for the 1500A emissivity. 

In the Bower06 model, ionizing emissivities are ~ 7 — 40 
times lower than in the Baugh05 model over the redshift range 
2 = 3 — 10 (assumi ng a 100% escape fraction in both cases). 
iRaicevic et alj d2010al) have made a more detailed comparison of 
the emissivities in these t wo models. 

IBenson et alj d2006h investigated reionization in the Baugh05 
and other GALFORM models using an analytical Stromgren sphere 
model for the growth of ionized regions in the IGM, and found that 
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Figure 17. Rest-frame 1500A luminosity density (in comoving units) vs redshift in the Baugh05 model, (a) The blue lines show the far-UV luminosity 
density due to all galaxies, while the green and red curves show the separate contributions of quiescent galaxies and ongoing bursts. Solid and dashed curves 
show luminosity densities with and without dust extinction, (b) The contribution to the dust-extincted far-UV luminosity density from galaxies with 1500A 
luminosities brighter than different limits, as indicated in the key. 
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Figure 18. Ionizing emissivity (number of Lyc photons per unit time per unit comoving volume) vs redshift in the Baugh05 model, (a) The blue curve shows 
the total emissivity (without dust extinction) and the green and red curves show the separate contributions of quiescent galaxies and ongoing bursts, (b) The 
contribution to the ionizing emissivity from galaxies with 1500A luminosities brighter than different limits, as indicated in the key. 



reionization was predicted to occur at z ~ 12 for a Lyc escape frac- 
tion of 1 00%. Reioniza t ion in t his model is studied in much greater 
detail in iRaicevic et aT] ( l2010j|bl) . who use a radiative transfer sim- 
ulation to show that, for a Lyc escape fraction in the range 0.1-1, 
reionization should occur at z ~ 8 — 10. 



7 CONCLUSIONS 

We have made a detailed investigation of the properties and evo- 
lution of Lyman-break galaxies (LBGs) predicted by hierarchi- 
cal models of galaxy formation. We followed the galaxy forma- 
tion process in the framework of the ACDM cosmology using the 
GALFORM semi-analytical model, which includes physical treat- 
ments of the hierarchical assembly of dark matter halos, shock- 



heating and cooling of gas, star formation, feedback from super- 
nova explosions, AGN and photoionization of the IGM, galaxy 
mergers and chemical enrichment. The luminosities of galaxies are 
calculated from a stellar population synthesis model, and dust ex- 
tinction is then included using a self-consistent theoretical model 
based on the results of radiative transfer calculations. The dust mass 
is calculated from the predicted mass and metallicity of the cold gas 
component, and this is combined with the predicted galaxy radius 
to calculate the dust extinction optical depth. The far-UV dust ex- 
tinction is a critical component in any model for LBGs. 

We have presente d prediction s for two variants of the 
GALFORM model. In the lBaugh et alj (2005, Baugh05) model, the 
formation of very massive galaxies is inhibited by supemova-driven 
superwinds which eject gas from halos, star formation at high red- 
shifts is dominated by starbursts triggered by galaxy mergers, and 
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stars form in these bursts with a top-heavy IMF. The top-heavy IMF 
was motivated by the need to explain the number counts and red- 
shift distributions of the faint sub-mm galaxies. This model also 
matches a wide range of other d ata on local galaxies (such as 
gas masses and disk sizes). In the iBower et alj (2006, Bower06) 
model, the formation of very massive galaxies is instead inhibited 
by AGN feedback which heats the gas in halos, starbursts play a 
much smaller role in star formation, and all stars form with a Solar 
neighbourhood IMF. The Bower06 model underpredicts the sub- 
mm number counts by more than an order of magnitude, due to 
having too few very luminous and dusty star-forming galaxies at 
high redshifts. This shortcoming might be remedied by introduc- 
ing a top-heavy IMF in bursts, as we will explore in a future paper. 
Both models match the present-day optical and near-IR luminosity 
functions. 

We first considered the evolution of the rest-frame far-UV 
luminosity function (f|3j. Both models predict modest evolution 
over the redshift range 3 < 2 < 8, but a more rapid evolution 
at higher redshifts, z > 8, driven by the build-up of dark mat- 
ter halos. However, the models differ in that the Bower06 model 
predicts a more extended high-luminosity tail than the Baugh05 
model, once dust extinction is included. The effects of dust ex- 
tinction on the far-UV luminosity function are much larger in the 
Baugh05 model (~ 2 mag) because the bright end of the lumi- 
nosity function is dominated by starbursts in which the dust con- 
tent is enhanced by metal production with the top-heavy IMF. We 
made a detailed comparison of the predictions of both models with 
observed far-UV luminosity functions of LBGs over the redshift 
range z — 3 — 10. We found that the Baugh05 model, without any 
modification of its parameters, predicts a far-UV luminosity func- 
tion in excellent agreement with current observational data over 
this whole range. On the other hand, the Bower06 model conflicts 
with the LBG observations at z = 3 — 7 because it predicts too 
many high-luminosity galaxies. We then investigated the effect on 
the predicted luminosity functions of varying some of the model 
parameters from their default values. Assuming a Solar neighbour- 
hood, rather than top-heavy, IMF in bursts has only a modest im- 
pact on the far-UV luminosity function, because the effects of lower 
intrinsic stellar luminosities are partly compensated by lower dust 
extinctions. However, such a model predicts far too few sub-mm 
galaxies. We find that the luminosity function over the range cov- 
ered by observational data is fairly sensitive to the assumed star 
formation timescale in bursts, especially at higher redshifts, but is 
less sensitive to the strength of supernova feedback, when these pa- 
rameters are varied over physically reasonable ranges. The inability 
of the Bower06 model to match the observed LBG luminosity func- 
tion data appears to be caused mainly by the short assumed star for- 
mation timescale in bursts, rather than by the AGN feedback model 
or the assumed IMF. 

We next investigated a wide range of other physical proper- 
ties of LBGs predicted by the models (Q. We first considered the 
sizes of galaxies in the rest-frame far-UV, and compared to recent 
observational measurements at z ~ 2 — 7. We found that both 
models predicted sizes in reasonable agreement with observations 
at higher UV luminosities, but only the Baugh05 model is consis- 
tent with observed sizes at lower luminosities. We then presented 
predictions of the Baugh05 model for stellar, halo and gas masses, 
clustering bias, circular velocity, burst fractions, bulge-to-disk ra- 
tios, star formation rates and gas and stellar metallicities, and made 
brief comparisons with relevant observational data. The model pre- 
dictions appear to be broadly compatible with current observational 
constraints (many of which are rather uncertain) in most cases. A 



particularly interesting issue is the stellar masses - our predicted 
values are well below observational estimates based on fitting stel- 
lar population models to broad-band photometry. However, the ob- 
servational estimates all assume a Solar neighbourhood IMF, while 
in the Baugh05 model the LBG population is dominated by star- 
bursts forming stars with a top-heavy IMF. The observationally in- 
ferred stellar masses therefore cannot be directly compared with 
the values from the model. When instead we compare the IR fluxes 
(which drive the observational stellar mass estimates) directly, the 
model is much closer to the observations. We will investigate this 
important issue in more detail in a future paper. We will also make 
a more detailed study of LBG clustering in future work, since this 
provides constraints on the masses of the dark matter halos hosting 
LBGs. 

In §5\ we showed predictions for LBGs at very high redshifts 
(z = 7 — 20) in the Baugh05 model, including surface densities 
of objects down to very faint apparent magnitudes (rriAB = 32), 
relevant for observations with future telescopes such as JWST and 
ELTs. We find that deep surveys planned with JWST should be able 
to detect a few LBGs at z ~ 15 and ttiab ~ 30 — 31, and many 
more at lower redshifts. LBGs detected at ulab ~ 31 are predicted 
to have angular radii ~ 0.02 — 0.05 arcsec, depending only weakly 
on redshift over this range, and to have circular velocities ~ 30 — 
100km s -1 , again only weakly dependent on redshift. 

In fj6] we showed the predicted evolution of the far-UV lumi- 
nosity density, and its relation to the cosmic SFR history, again 
in the Baugh05 model. The unextincted 1500A luminosity den- 
sity tracks the SFR density in high-mass stars (m > 5Mq) more 
closely than the total SFR density, since the relative contributions 
of quiescent and burst star formation (with Solar neighbourhood 
and top-heavy IMFs respectively) change with redshift. The effect 
of dust extinction on the far-UV luminosity density is predicted to 
be large, ~ 2 mag at 1500A in the range 3 < z < 15, dropping 
to ~ 1 mag at z = 0. Finally, we considered the predicted contri- 
bution of galaxies to the emissivity of ionizing photons which can 
reionize the IGM. For a constant escape fraction of ionizing pho- 
tons from galaxies, this emissivity falls by a factor ~ 100 from its 
peak at 2 ~ 5 to 2 = 15. At high redshift, most of the ionizing 
photons are predicted to come from very low luminosity galaxies, 
so that, for example, to detect the galaxies responsible for > 50% 
of the ionizing emissivity at 2 = 10 would require an LBG sur- 
vey probing fainter than A/ab(1500A) — 5 log h ~ —15, corre- 
sponding to an apparent magnitude ttiab ~ 32 at the same rest- 
frame wavelength. The predictions of our model f or reionization 
of the I GM are discussed in much greater detail in lRaicevic et al.l 
teOlOalbl) . 

In conclusion, we find that the iBaugh et al.l d2005l) model, 
which was originally constructed to match the far-UV luminos- 
ity function of LBGs only at 2 = 3, predicts results in remark- 
ably good agreement with subsequent observations of LBGs out to 
2 = 10. Further exploration of whether this model provides a phys- 
ically accurate description of LBGs and other high-redshift galaxy 
populations will require more detailed comparisons between the 
model predictions and observational data (e.g. for stellar masses, 
clustering, and colours), but also new and more sensitive observa- 
tions. 
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APPENDIX A: COMPARISON OF THE BAUGH ET AL. 
AND BOWER ET AL. MODELS 

The Baugh05 and Bower06 models were originally constructed 
with somewhat different aims, and their parameters for star for- 
mation, feedback etc were chosen by comparison with differ- 
ent, though partly overlapping, observational datasets. The pa- 
rameters for both models were constrained to reproduce the ob- 
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served present-day B- and K-band galaxy luminosity functions. 
The Baugh05 model was constructed to also match a wide range 
of other properties of present-day galaxies, including the far-IR 
(60/im) luminosity function and the dependence of disk size, gas 
content and metallicity on galaxy luminosity. The Bower06 model 
did not use any of these present-day properties as constraints in 
choosing its parameters, but instead aimed to reproduce the gen- 
eral form of the colour-magnitude distribution of galaxies at z — 0. 
In addition, the Baugh05 model was designed to explain the ob- 
served number counts and redshift distributions of the faint sub-mm 
(850/im) galaxies, and the rest-frame far-UV luminosity function 
of Lyman-break galaxies at z ~ 3, while the Bower06 model in- 
stead focused on matching the observed evolution of the K-band 
luminosity function of galaxies at z < 2. 

Given these differences, it is not surprising that the two mod- 
els, although fundamentally quite similar, have dif ferent successes 
and dr awbacks, as discussed in subsequent papers. Gonzal ez et al.l 
J2009h compared both models side-by-side with observational data 
from the SDSS on luminosity functions, colours, morphologies and 
sizes of present-day galaxies. They found that both models are 
reasonably consistent with the fraction of early-type galaxies as 
a function of luminosity, but that the Bower06 model is in bet- 
ter agreement with the distribution of galaxy colours. Further, the 
Baugh05 model agrees much better with the size-luminosity rela- 
tion of disk galaxies, but both models have proble ms explaining 
the siz e s of spheroidal galax ies (as also found by lAlmeida et al.l 
d2007h ). lAlmeida et al] J2008h found the Bower06 model in better 
agreement with the luminosity function and co rrelation function of 
Luminous Red Galaxies (LRGs) in the SDSS. IPower et al.l j201Ch 
investigated the cold gas mass functions in both models, and found 
the Baugh05 model in better agreement with data from the HIPASS 
21cm survey of the local universe. 

At higher redshifts, the Baugh05 model seems in generally 
better agreement with observations of actively star-forming galax- 
ies at z > 2, while the Bower06 model is in bette r agreement 
with o bservations of more passive galaxies at z < 2. iBower et al.l 
J2006I) compared the evolution of the K-band luminosity function 
in both models with observations at z ^ 1.5, and found that 
the Baugh05 model predicts somewhat too few high-luminosity 
galaxies at higher redshifts. The Baugh05 model also predicts 
much lower number counts for Extremely Red Objects (EROs) 
than Bower06, with th e latter being in much be tter agreement 
with observational data dGonzalez-Perez et alj2009l) . However, the 
Bower06 model predicts number counts of faint sub-mm galaxies 
which are too low by more than an order of magnitude. The com- 
parison of both models with observations of LBGs is the topic of 
this paper. 



© 0000 RAS, MNRAS 000, 000-000 



